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.1 

Kinetics  of  crack  propagation  in  hydrogen  charged  niobium  were 

studied  in  the  temperature  range  77  to  303  K  and  for  hydrogen 

concentrations  between  0.068  and  3.2  at.%.  Crack  velocities  between 
-8  -2 

10  and  10  m/s  were  determined  as  a  function  of  the  Mode  I  stress 
intensity  factors  using  single  edge  notch  and  double  torsion  specimens. 
Validity  of  the  double  torsion  technique  was  established  using  a  variety 
of  methods  including  photoelastic  analysis  of  the  stresses  in  the  plane 
normal  to  the  specimen's  crack  surface. 

Increasing,  decreasing  and  constant  K^.  tests  were  carried  out  and 
the  InV  versus  k_j.  data  exhibited  various  behavior  depending  upon  the 
type  of  loading.  Linking  up  of  fracturing  hydrides  was  shown  to  be 
responsible  for  a  unique  three  region  (V-K)  behavior  which  occurred 
only  at  intermediate  temperatures  for  decreasing  K^  tests.  Analysis 
of  the  temperature  dependence  for  Stage  II  during  constant  tests 
indicated  that  the  cracking  kinetics  are  controlled  by  hydrogen  diffu¬ 
sion  to  the  crack  tip.  Differences  between  crack  velocities  in  deuterium 
and  hydrogen  charged  specimens  were  consistent  with  the  isotopic 
differences  in  their  diffusion  coefficients.  Fracture  surfaces  of  the 
test  specimens  were  examined  and  shown  to  exhibit  features  consistent 
with  their  crack  velocity  behavior. 
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1.  INTRODUCTION 


The  phenomenon  of  hydrogen  embrittlement  of  pure  metals 
and  alloys  has  been  extensively  researched. ^  It  is  now 
recognized  that  a  single  mechanism  cannot  account  for  all  instances 
of  hydrogen  embrittlement.  Based  on  the  degree  of  hydrogen  solubility 
in  metals,  two  classes  of  hydrogen  embrittlement  can  be  distinguished. 
The  first  class  includes  metals  that  have  low  hydrogen  solubility, 
endothermic  heats  of  solution,  and  do  not  form  stable  hydrides, 
e.g.  copper,  iron  and  molybendum.  Some  elements  in  this  class, 
e.g.  iron,  nickel  and  their  alloys  have  severe  susceptibility  to 
hydrogen  embrittlement,  while  others  e.g.,  aluminum  and  molybendum 
have  moderate  susceptibility  or  are  insensitive  like  copper.  The 
second  class  includes  metals  that  have  high  hydrogen  solubility, 
exothermic  heats  of  solution  and  form  stable  hydrides.  This  class 
includes  all  elements  belonging  to  groups  IVa  and  Va  of  the  periodic 
table  and  many  others  which  form  stable  hydrides  such  as  lithium, 
titanium  and  magnesium.  The  formation  of  stable  hydrides  in  these 
metals  causes  their  severe  susceptibility  to  hydrogen  embrittlement. 
Niobium  is  a  member  of  the  second  class. 

In  specimens  which  form  stable  hydrides  while  under  external 
stress,  the  embrittlement  process  involves  two  main  steps.  The  first 
step  requires  hydrogen  to  diffuse  down  a  stress  gradient  to  the 
stress  concentrator  where  a  hydride  is  formed.  The  second  step  is 
the  fracture  of  the  hydride.  The  crack  tip  advances  by  a  repetition 
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of  these  steps.  The  first  step  controls  the  crack  velocity  since 
the  hydride  itself  fractures  quite  rapidly. 

This  study  investigates  the  kinetics  of  crack  propagation 
in  hydrogen  embrittled  niobium  to  further  understand  the  embrittlement 
mechanism.  Kinetics  are  presented  in  terms  of  crack  velocity  as  a 
function  of  stress  intensity  factor.  The  fracture  mechanics  approach 
was  used  since  the  stress  intensity  can  relate  the  crack  tip  stresses 
to  the  externally  applied  stresses.  Double  torsion  and  single  edge 
notch  specimen  were  used  in  crack  velocity  measurements. 

Photoelastic  analysis  of  the  double  torsion  specimen  was  undertaken 
and  evidence  for  Mode  I  stress  intensity  at  the  crack  tip  is 
presented.  Fractography  was  used  to  correlate  fracture  surface 
features  to  crack  velocity  kinetics.  Crack  velocity  measurements  as  a 
function  of  increasing,  decreasing  and  constant  stress  intensity 
were  made.  Several  contributing  factors  are  shown  to  explain  the 
crack  velocity  behavior  in  each  case.  These  factors  include  the 
type  of  loading  a  specimen  undergoes,  the  magnitude  of  the  loading 
and  the  test  temperature  proximity  to  the  specimen's  solvus 
temperature.  The  ideas  presented  here  are  generally  applicable  to 
other  exothermic  hydrogen  occluders  forming  stable  hydrides. 


i. 


2.  REVIEW  OF  PREVIOUS  RESEARCH 

2.1  Stress  Effects  on  Hydrogen  Solubility 

An  expression  relating  the  solubility  of  hydrogen  in  a  stress 

free  solid  to  the  solubility  in  a  stressed  solid  was  given  by 

(2) 

Grossbeck  as 

C°/c£  =  exp(WH  -  wH)/RT  (1) 

where  W  =  total  work  per  mole  addition  of  hydrogen, 

H 

w„  *  partial  molar  strain  energy, 
n 

Equation  (1)  was  derived  using  the  thermodynamic  approach  of  stressed 

solids  of  Li  et  al.  ^  Puls^  has  shown  an  identical  result  can 

be  derived  for  hydride  forming  alloys  using  the  atomistic  model  of 

(2  5) 

Eshelby.  It  was  further  shown  ’  that  for  a  specimen  subjected 
to  a  triaxial  tension  stress  field  and  assuming  all  shear  stresses  to 
be  zero,  Eqn.  (1)  becomes 

CK  ■  “p'V  VRT)  <2) 

where  V  -  partial  molal  volume  of  hydrogen  in  the  metal, 

]?hyd  *  hydrostatic  stress  given  by  (a.^  +  a22  +  o33)/ 3. 

Values  calculated  for  C„/C°  in  zirconium  ranged  from  1.65  to  1.005 

ri  H 

depending  upon  temperature  and  stress. ^  For  niobium,  in  plane 
stressed  specimens,  Grossbeck^^  has  shwon  the  ratio  /CS  co  be 
as  high  as  66  for  a  temperature  of  77  K  and  a  stress  of  one  tenth 


the  shear  modulus. 
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2.2  Stress  Effects  on  Hydride  Solvus 

(2) 

Grossbeck  noted  that  the  increase  of  hydrogen  at  a  crack 
tip  from  an  applied  stress  does  not  by  itself  result  in  hydride  pre¬ 
cipitation  for  an  exothermic  occluder.  The  chemical  potential  of  the 
hydride  phase  must  be  lowered  upon  application  of  an  external  stress  in 
order  for  the  solvus  to  be  raised.  Flanagan  et  al. ^  have  dis¬ 
tinguished  between  the  effects  of  uniform  and  nonuniform  applied  stress 
fields  on  the  thermodynamics  of  hydride  precipitation.  Elastic  and 
plastic  constraints  were  not  considered.  Their  conclusions  can  be 
summarized  as  follows. 

1.  For  systems  where  partial  molal  volumes  of  hydrogen 
in  the  hydride  and  solid  solution  are  equal,  a 
uniform  stress  has  no  effect  on  the  solvus. 

2.  For  the  case  where  the  partial  molal  volumes  are  not 
equal,  a  uniform  stress  has  a  small  effect  on  the 
solvus  depending  upon  the  magnitude  to  the  differences 
in  partial  molal  volumes. 

3.  For  a  nonuniform  stress,  a  hydride  will  precipitate  when 
the  chemical  potentials  of  the  unstressed  volume  of 

the  solid  solution  and  the  stressed  volume  of  the 
hydride  are  equal.  This  occurs  even  in  the  absence 
of  a  structural  change  in  forming  the  hydride. 

Westlake  and  Ockers,^  using  electrical  resistivity  methods, 
have  shown  that  the  terminal  solid  solubility  for  the  niobium-hydrogen 
system  depended  on  the  cooling  rate.  In  addition,  the  alpha-beta 
solvus  exhibited  a  hysteresis  for  the  precipitation  and  reversion 
of  the  beta  phase.  The  magnitude  of  the  hysteresis  depended  on  the 
number  of  thermal  cycles  the  specimen  underwent  and  in  the  absence 
of  an  external  stress,  they  found  the  solvus  to  be  given  by 
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In  C°  =  1.320  -  11.45(kJ/mol)/RT°(n)  (3) 


where  C°  »  hydrogen  concentration  in  atom  fraction, 
rl 

T°(n)  *  stress  free  solvus  temperature  after  many  thermal  cycles. 

It  has  long  been  known  that  the  formation  of  niobium  hydride 

/Q\ 

is  accompanied  by  a  12%  increase  of  the  molar  volume.  This  led 
(9) 

Birnbaum  et  al.  to  consider  the  effects  of  elastic  and  plastic 
deformations  resulting  from  hydride  precipitation  in  the  alpha-beta 
solvus.  They  showed  that  the  elastic  accommodation  due  to  cooling 
was  reversed  during  re-solution  and  the  plastic  accommodation  was 
responsible  for  the  hysterisis  in  the  heating  and  cooling  cycles. 

Electron  microscopy  has  been  used  to  observe  the  accommodation 
effects  of  hydride  precipitation  for  niobium-nvdrogen  alloys.  ^ 

It  was  found  that  during  precipitation,  volume  increase  of  the  beta  phase 
was  accommodated  elastically  and  by  punching  out  cf  interstitial 
dislocation  loops  in  the  alpha  phase.  Upon  disolction  of  the  hydride 
during  heating,  volume  decrease  was  accommodated  nv  recovery  of 
elastic  strain  and  dislocation  processes. 

Grossbeck  and  Birnbaun/'^  showed  that  the  alpha-beta  solvus 
increases  in  the  presence  of  an  external  tensile  stress  and  its  free 
energy  change  is  given  by 


-  £0  . e  .  •  V 
L  ij  ij 


where 


applied  stress, 

strain  tensor  describing  the  alphi-aeti  chase  change. 


V  *  molar  volume. 
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It  was  also  shown,  that  a  further  increase  in  the  solvus  results  from 

the  external  stress  assisting  the  plastic  accommodation  process 

(2) 

during  precipitation.  Grossbeck  showed  that  the  solvus  temperature 
given  in  Eqn.  (3),  is  increased  in  the  presence  of  an  external 
stress  by  the  amount  given  by 

AT  =  -PAV/R(1.320  -  In  Cu)  (5) 

S  rl 

where  P  =  applied  hydrostatic  stress, 

AV  *  volume  change  accompanying  hydride  formation. 

For  applied  stress  of  0.05  of  the  shear  modulus  under  plane  stress 

conditions,  increases  of  18  K  to  28  K  were  found  for  hydrogen 

concentrations  between  0.1  and  2.0  at.  %,  respectively.  Birnbaum 
(9) 

et  al.  have  calculated  the  equilibrium  solvus  based  on  a 
consideration  of  the  plastic  and  elastic  contributions  to  the  molar 
free  energy  change  in  forming  the  beta  phase  hydride.  Puls,^ 
however,  contends  tha t  the  equilibrium  solvus  cannot  in  general  be 

(121 

calculated  but  must  be  measured  experimentally.  Koike  and  Suzuki 
presented  evidence  for  triaxial  stresses  n'u  contributing  to  stress 
induced  hydride  formation  at  a  crack  tip  in  the  case  of  vanadium-hydrogen 
alloys.  For  these  alloys,  the  alpha  to  beta  phase  transformation 
resembles  a  martensitic  transformation  and  a  uniaxial  tension  appears 
responsible  for  the  stress  induced  hydride  formation. 

2. 3  Subcritical  Crack  Growth  Behavior  and  Mechanisms 

It  has  been  observed^^  that  subcritical  crack  growth 
often  exhibits  three  distinct  stages  of  crack  velocity  (V)  versus 
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stress  intensity  (K^)  ,  (Fig.  1).  Stage  I  occurs  at  the  low  values 
and  low  crack  velocities  and  Ls  characterized  by  a  steep  dependence 

of  V  on  and  a  possible  threshold  stress  intensity,  .  Stage  II 

o 

is  characterized  by  a  plateau  region  where  the  crack  velocity  is 


independent  of  K^..  Stage  III  represents  the  onset  of  unstable  crack 


growth  and  rapid  fracture  near  K^. 


Several  investigations  ^  of  st 


ress  corrosion  cracxing 


exhibited  three  stage  (V-K)  curves.  It  was  suggested  that  Stage  II 

results  from  a  rate  limiting  mechanism  such  as  diffusion  of  the 

corrosive  component  to  the  crack  tip.  Three  stage  behavior  has  also 

(21) 

been  found  by  Nelson  for  titanium  alloys  fractured  in  hydrogen 
gas  atmospheres.  He  presented  a  film  rupture  model  of  hydrogen- 
induced  slow  growth  for  these  alloys  which  include  6A1-4V,  5Al-2.5Sn 
and  8AI-IM0-IV.  They  share  common  microstructures  all  of  which 
exhibit  acicular  alpha  phase  platelets  in  a  beta  phase  matrix. 

Three  stage  behavior  has  been  observed  for  certain  high  strength  stee 


e.g.  4130  and  4340  in  gaseous  hydrogen. 


Non-metallic  svstems 


have  also  exhibited  three  stages  of  crack  velocity.  Examples  of 

(23) 

these  include  polycrystalline  alumina  in  toluene, 

polymethylmethacrylate  in  air^“^  and  soda  lime  glass  in  air.^^ 

On  the  basis  of  the  behavior  of  18-Ni-maraging  steels  in 

f  1 8 ) 

gaseous  hydrogen,  Hudak  and  Wei  has  suggested  that  the 
independence  of  Stage  II  velocities  is  due  to  a  rate  limited  chemical 
process  which  restricts  the  amount  of  hydrogen  delivered  to  the 
fracture  zone.  An  activation  energy  for  the  rate  limiting  step 


can  be  calculated  for  constant  tests  using  a  family  of  Stage  II 


9 


velocity  data  taken  at  various  temperatures.  Activation  energies 
measured  for  these  steels  do  not  as  yet  suggest  a  definite  mechanism. 

In  particular  it  is  noted  that  activation  energies  for  hydrogen 
permeation  in  these  steels  are  two  to  three  times  larger  than  the 

( 18 ) 

activation  energies  which  characterize  Stage  II  crack  propagation. 

(24) 

Williams  and  Nelson  suggested  that  the  adsorption  of  hydrogen  on 

steel  surfaces  is  the  rate  limiting  step  of  gaseous  hydrogen  assisted 

crack  growth.  A  similar  rate  limiting  process  due  to  surface  transfer 

(25) 

is  advocated  by  Oriani  and  Josephic  in  their  lattice  decohesion 

model  for  hydrogen  assisted  cracking  in  steels. 

The  occurrence  of  three  stage  behavior  in  hydrogen  embrittled 

refractory  alloys  has  been  well  established.  Early  work  on  Zircalloy 
(26 ) 

by  Dutton  et  al.  '  suggested  Stage  I  and  Stage  II  were  independent 

phenomenon  being  controlled  by  elastic  and  plastic  zone  stress  fields, 

(27) 

respectively.  However,  more  recent  work  by  Simpson  and  Puls 
suggested  that  Stage  I  is  not  due  to  a  different  mechanism  than  Stage  II. 
Indeed,  Stage  II  does  not  even  appear  to  exist.  Also  much  emphasis 
was  placed  on  the  conditioning  of  the  machined  specimen  notch  in 
order  to  obtain  what  they  considered  valid  test  results.  It  will  be 
shown  in  this  thesis,  that  three  stage  behavior  for  hydrogen  embrittled 
niobium  exists  and  its  absence  or  presence  is  a  function  of  the  type 
of  loading  a  specimen  undergoes. 


m 
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2.4  Crack  Propagation  Measurements  and  Techniques 

Techniques  for  measuring  crack  propagation  fall  in  two 
categories.  The  first  requires  the  direct  visual  observation  of 
the  crack  tip  either  by  eye  or  video  tape  recorder.  The  second 
involves  indirect  measurement  of  the  crack  velocity.  This  would 
include  acoustic  emission,  potential  drop,  load  relaxation 
and  crosshead  displacement.  Each  of  these  techniques  have  advantages 
and  disadvantages.  The  type  of  equipment  available  to  an  investigator 
would  be  important  in  determining  which  technique  one  would  use. 

The  most  common  fracture  mechanics  specimens  used  in  crack 
propagation  measurements  are  the  single  edge  notched  (s.e.n.), 
double  cantilever  beam  (d.c.b.),  center  cracked  tensile  (c.c.t.), 
and  compact  tension  (c.t.)  specimens.  Less  used  but  gaining  in 
popularity,  for  reasons  explained  below,  is  the  double  torsion 
specimen  (d.t.). 

Only  the  tapered  d.c.b.  specimen  and  d.t.  specimen  are 
constant  specimens.  That  is,  is  independent  of  the  crack  length. 
This  is  particularly  useful  since  K^.  is  then  simply  proportional  to 
the  applied  load.  The  measurements  of  crack  velocities  in  specimens 
other  than  the  d.t.  specimen  are  more  difficult.  The  potential  drop 
method  when  applied  to  the  c.t.,  tapered  d.c.b.  or  wedge  opening 
loaded  etc.  specimens  requires  very  stable  and  low  noise  constant 
current  supplies.  Careful  calibration  techniques  are  required  when  a 
change  in  load  is  attempted  during  a  test  In  order  to  relate  the 
potential  drop  to  the  correct  crack  extension. 


Crack  velocity  measurements  in  the  d.t.  specimen  are 
particularly  simple.  As  described  in  the  Theory  section  in  more 
detail,  only  the  monitoring  of  either  ’"he  load  relaxation  or 
crosshead  displacement  as  a  function  of  time  is  required.  Moreover, 
the  geometry  of  the  d.t.  specimen  consists  simply  of  a  thin  plate, 
minimizing  the  amount  of  machining  and  metal  required  for  a  specimen. 
It  was  primarily  for  these  reasons  that  the  d.t.  specimen  was  chosen 
for  these  test. 

The  d.t.  specimen,  whose  geometry  is  shown  in  Fig.  2  was 

(29) 

designed  by  Outwater  and  Gerry  but  first  reported  by  Kies  and 
Clark.  It  was  developed  in  response  to  the  need  for  a  specimen 

which  was  simple  to  machine  and  enable  crack  velocities  to  be 
measured  in  severe  testing  environments.  The  equations  governing 
the  crack  velocity  and  stress  intensity  for  the  d.t.  specimen  are 
derived  in  the  Theory  section.  These  equations,  however,  apply 
strictly  to  specimens  undergoing  elastic  deflection  only.  This 
assumption  is  true  for  glass  or  ceramics  (at  room  temperature)  but  is 
not  generally  correct  for  steels  or  niobium.  Another  assumption 
in  the  derivation  of  these  analytic  equations  involves  the  crack  front 
profile.  Ideally,  the  crack  front  is  normal  to  the  specimen's  top 
surface  whereas  in  reality  it  is  a  sloping  curve  as  shown  in 
Fig.  3.  Departure  from  a  perpendicular  crack  front  increases  as 
one  goes  to  increasingly  ductile  materials.  To  correct  for  the 
nonideal  crack  front  profile,  a  multiplicative  constant  derived 
in  section  5.2.1  is  applied  to  the  analytical  equations.  This  factor 
must  be  empirically  found  for  each  material  but  does  not  differ 
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Fi^.  3.  Tvpical  crack  front  profile  of  d.t.  specimen. 
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(31) 

greatly  from  unity  in  most  cases.  Pollet  and  Bums  and  Shettv 
(32) 

et  al.  have  derived  additional  crack  velocity  correction  factors 

for  the  d.t.  specimen  which  essentially  give  the  same  results  as  the 

method  used  in  this  investigation. 

(33) 

Trancina  measured  the  stress  intensity  factor  of  the  d.t. 

specimen  as  a  function  of  crack  length  and  ligament  length  using  a 

three  dimensional  elastic  finite-element  stress  analysis.  He  found 

the  stress  intensity  was  nearly  constant  through  one  half  the 

specimen's  thickness  and  the  analytical  equations  were  correct  for 

crack  lengths  greater  than  0.55  W  and  ligament  lengths  less  than  0.65W 

where  W  is  the  specimen's  width. 

An  additional  test  for  the  validity  of  the  d.t.  specimen 

is  to  compare  its  results  with  measurements  taken  from  specimen 

geometries  which  are  known  to  give  "correct"  fracture  mechanics  data. 

A  comparison  between  bend  specimens  and  d.t.  specimens  agreed  very 

(23) 

well  in  measurements  of  the  critical  stress  intensity  factor. 

Comparisons  of  the  (V-K)  data  for  glass  and  4130  steel  using 

d.c.b.  and  d.t.  specimens,  showed  very  similar  functional  dependence 

except  that  the  d.t.  specimens  had  a  consistently  higher  crack  velocity 

(23) 

by  a  factor  of  five.  This  discrepancy  has  been  attributed  to 

the  nonideal  crack  front  profile  discussed  above.  For  niobium, 
during  constant  tests,  the  bottom  crack  advance  length  measured 
after  a  run  with  a  micrometer,  was  within  10%  of  the  length  calculated 
from  velocity  and  time  data.  Various  other  checks  of  the  validitv 
of  the  d.t.  specimen  geometry  fo  fracture  of  the  niobium-hydrogen 


alloys  were  made  as  will  be  discussed  and  these  all  indicated  that 
the  tests  give  valid  data.  Procedures  for  taking  into  account 
the  contribution  to  the  deflection  from  the  plastic  behavior  of 
niobium  will  be  presented  in  the  Experimental  Procedure  section. 
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3.  THEORY 

The  usefulness  of  the  d.t.  specimen  geometry  in  measuring 
subcritical  crack  growth  as  a  function  of  stress  intensity  factor 
is  due  to  two  properties.  One  is  that  the  stress  intensity  is 
independent  of  the  crack  length,  i.e.  a  constant  specimen. 
Secondly,  the  crack  velocity  can  be  calculated  by  monitoring  either 
the  rate  of  load  relaxation  (constant  crosshead  deflection  test)  or 
the  rate  of  crosshead  deflection  (constant  load  test) .  The  deriva¬ 
tion  of  the  equations  governing  the  velocity  at  constant  load,  V  , 
and  at  constant  crosshead  deflection,  V^,  is  given  below. 

3 . 1  Constant  Deflection  Double  Torsion  Test  (Decreasing  K^.) 

(34) 

Novozhilov  showed  that  the  d.t.  specimen's  compliance, 
C,  for  small  deflections  and  for  plates  whose  width,  W,  is  much 
greater  than  its  thickness,  T,  is  given  by 

C  =  (3w  a/WT3u)  =  y/P  (6) 

m 

where  a  »  crack  length  measured  from  specimen’s  edge. 

w  »  moment  arm  length, 
in 

U  »  shear  modulus, 

P  »  load, 

y  ■  deflection  of  specimen  due  to  loading. 

(23) 

Williams  and  Evans  showed  that  the  stress  intensity  factor. 


n 

assuming  the  crack  front  profile  is  independent  of  crack  length, 
is  given  by 

Kt  *=  Pv  [3(1  +  v)/WT3T  ]L/2  (7) 

I  m  n 

where  v  *  Poisson's  ratio, 

T  «  Thickness  of  specimen  minus  groove  depth, 
n 

Since  the  factors  v,  W,  T  and  T  are  constants  which  depend  only  on 

n 

the  specimen's  dimensions  and  mechanical  properties,  they  can  be 
combined  into  one  constant  to  yield 

C,  *  W  [3(1  +  v)/WT3T  ]1/2  (8) 

l  m  n 

-*3/2 

The  constant,  C^,  has  units  of  (length)  .  Equation  (7)  can  then 
be  rewritten  as 

K  =  C1P  .  (7) 

For  constant  deflection  (decreasing  K^) ,  a  =  a(t),  P  =  P(t)  and 
y  *  constant.  Rearranging  Eqn.  (6),  there  follows 

y (t)  «  3w2aP/WT3U  (6') 

m 

and  since 

U  -  E/2 (1  +  v)  (9) 

where  E,  is  Young's  modulus,  implies 

y (t)  -  6w2(1  +  v)aP/WT3  E  . 


(10) 
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Differentiating  Eqn.  (10)  with  respect  to  time,  t,  yields 

Oy/9t)  -  0  -  [3w2/WT3y]  •  [  (3a/3t)  P  +  (3P/3t)  a]  (11) 

y  m  y  y 

which  implies 

(3a/3t)y  =  -(a/P)  •  OP/St)^  .  (12) 

From  Eqn.  (6) , 

a  =  WT3Ey/6Pw2(l  +  v)  (13) 

ra 

and  substituting  the  above  expression  for,  a,  into  Eqn.  (12)  yields 

V  =  (3a/3t)  =  -(WT3Ev/P26w2(l  +  v) ]  •  (3P/3t)  .  (14) 

y  y  m  y 

Again  grouping  the  constants,  there  follows, 

C.  =  -WT3E/6w2(l  +  v)  (13) 

i.  m 

where  C2  has  units  of  force.  Equation  (14)  can  then  be  rewritten  as 
V  =  (C.,y/P2)  •  (3P/3t)  .  (14') 

y  2  y 

3 . 2  Constant  Load  Double  Torsion  Test  (Constant  K^) 

For  constant  load  conditions,  y  =  y(t),  a  *=  a(t)  and  P  = 
constant.  Letting 

a  -  6w2  (1  +  v)/WT3E  ,  (16) 

ra 


Eqn.  (10)  becomes 
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Differentiating  the  above  expression  with  respect  to  tint 

(3a/3t)  -  (aP)_1(3y/3t) 

P  P 

or 

V  =  (3a/3t)  =  [WT3E/6w2(1  +  v)P)  •  (3y/3t)  . 

ppm  p 

Expressing  the  above  equation  in  terms  of  C2  results  in 
V  -  -(CJ P)  •  (3y/3t)  . 

p  z  p 

The  crack  length,  a,  is  taken  as  the  longest  por:  :  r. 
occurs  on  the  bottom  of  the  specimen.  The  consequence  01  tr. 
shown  in  the  compliance  measurements  of  the  d.t.  tests  in  tr. 
Results  Section.  Photoelastic  analysis  of  the  d.t.  geonetrv 
presented  in  the  Appendix.  It  will  be  shown  that  the  stress 
intensity  factor  given  in  Eqn.  (7),  closely  resembles  the  st 
intensity  for  Mode  I  loading,  i.e.  K^. 


:  s  j 


. 
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4.  EXPERIMENTAL  PROCEDURE 

4.1  Sample  Preparation 

4.1.1  Single  Edge  Notch  Specimens 

Niobium  specimens  having  a  single  edge  notch  geometry  were 
prepared  from  0.1  cm  rolled  sheet  obtained  from  Wah  Chang  Albany 
Corporation.  Typical  specimen  dimensions  are  shown  in  Fig.  4. 
Reference  lines  from  150  to  90  microns  apart  were  engraved  on  the 
specimen's  surfaces  with  a  diamond  stylus  (Fig.  5).  These  lines 
allowed  the  average  velocity  of  the  crack  to  be  calculated  by  the 
following  expression, 


V  =  Ax/ At  (19) 

where  Ax  *  crack  advance  distance. 

At  =  elapsed  time  of  crack  advance. 

To  remove  surface  contamination,  the  specimens  were  etched  after 

machining,  in  a  solution  with  equal  parts  by  volume  of  HNO^,  HF  and 

lactic  acid.  The  mass  spectrographic  analysis  before  annealing  is 

shown  in  Table  1.  All  s.e.n.  specimens  were  annealed  prior  to 

—8 

hydrogen  charging  in  a  vacuum  of  6.5x10  Pa  and  held  at  2300  K  for 
several  hours.  The  typical  grain  size  diameter  after  annealing 
was  0.3  cm.  Specimens  were  charged  with  hydrogen  concentrations 
from  0.068  to  0.74  at  Z  in  an  ultra  high  vacuum  charging  apparatus 
which  allowed  equilibration  of  the  specimens  with  hydrogen  gas 


TABLE  1 

Mass  Spectrographic  Analysis  of  Niobium 
(Concentrations  in  atomic  ppm) 


F 

S 

Cl 

Ti 

V 

Cr 

Mn 

Fe 

Cu 

Zn 

Zr 

Mo 

Ta 


0.4 

4 

3 

0.1 

0.2 

0.04 

0.03 

0.03 

0.1 

0.06 

0.06 

0.04 

10 


W 
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(purified  by  a  palladium  membrane  cell)  at  controlled  temperature 
and  pressure.  Vacuum  fusion  analysis  was  used  to  measure  the 
interstitial  impurities.  Typical  values  are  found  in  Table  2. 


A. 1.2  Double  Torsion  Specimens 

Double  torsion  specimens  were  prepared  from  0.2  cm  plate  obtained 
from  Wah  Chang  Albany  Corporation.  Mass  spectrographic  analyses  for 
the  substitutional  impurities  in  the  niobium  stock  are  shown  in 
Table  3.  The  ASTM  grain  size  was  4.5  in  the  rolling  direction. 

Electron  discharge  machining  was  used  to  cut  grooves  down  the 
center  of  one  side  to  a  depth  of  approximately  one  half  the  specimen 
thickness  and  was  also  used  to  cut  a  notch  to  a  length  of  3  cm  on  one 
end  of  each  specimen.  This  was  done  to  avoid  end  effects  which 
exhibit  nonlinear  behavior  of  the  compliance  change  as  a  function  of 
crack  length.  This  is  discussed  further  in  section  5.2.1.  Thirty 
gauge  wire  was  used  in  all  e.d.m.  Typical  dimensions  of  the 
specimens  after  machining  are  shown  in  Fig.  2.  Etching  of  the  d.t. 
specimens  were  carried  out  as  described  for  the  s.e.n.  specimens  after 

which  they  were  either  annealed  or  used  as  received.  Specimens  to 

-4 

be  annealed  were  placed  in  a  vacuum  of  10  Pa  and  held  at  2300  K  for 
several  hours.  Typical  final  grain  size  diameter  was  0.7  cm. 

The  o^arging  procedure  was  similar  to  the  one  described  tor 

s.e.n.  specimens  and  consists  of  heating  the  d.t.  specimens  in  ')  i  t  K 

-4 

under  a  vacuum  of  1.3  x  10  Pa.  Purified  hydrogen  gas  was  i  lien 
introduced  at  the  appropriate  pressure  to  achieve  the  desired  livtcyn 
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TABLE  2 

Vacuum  Fusion  Analysis  for 
Hydrogen  Charged  S.E.N.  Specimens 
(Concentrations  in  at.  %) 


Specimen 

0 

H 

N 

112 

0.064 

0.068 

0.028 

114 

0.035 

0.178 

0.019 

119 

0.090 

0.50 

0.016 

121 

0.080 

0.50 

0.023 

124-129 

0.076 

o 

Uo 

O 

0.025 
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TABLE  3 

Mass  Spectrogranhic  Analysis  for 
Niobium  Double  Torsion  Specimens 
(Concentrations  in  atomic  ppm) 


F 

10 

A1 

* 

3 

S 

8 

K 

3 

Ca 

15 

Cr 

2 

Fe 

2 

Ni 

3 

Cu 

1 

Ta 

30 

W 

10 

concentration  in  solution.  This  concentration  was  maintained  bv 
slowly  lowering  both  the  specimen  temperature  and  hydrogen  gas 
pressure  by  the  appropriate  amounts.  This  procedure  i ■;  followed 
until  a  temperature  of  approximately  400  K  was  reached.  It  is 
fortunate  that  an  oxide  layer  forms  on  the  surface  for  temperatures 
under  about  520  K  which  prevents  hydrogen  from  diffusing  out  after 
charging.  Deuterium  gas  was  used  to  charge  several  d.t.  specimens 
to  determine  whether  an  isotope  effect  could  be  seen  in  the  (V-K) 
curves.  The  charging  procedure  for  deuterium  was  the  same  as  for 
hydrogen  except  the  equilibration  times  are  longer  due  to  deuterium's 
smaller  diffusion  coefficient.  Vacuum  fusion  analysis  for  selected 
annealed  and  unannealed  specimens  are  shown  in  Table  4. 

After  charging,  the  d.t.  specimens  were  then  suitable  for 

testing  using  the  constant  deflection  method  (decreasing  K  ). 

Specimens  to  be  used  for  the  constant  load  tests  (constant  K  )  were 

fatigued  to  initiate  a  small  sharp  crack  in  front  of  the  machined 

notch.  Fatiguing  was  carried  out  on  a  four  column  table  MTS 

machine  using  a  frequency  of  2  Hz  and  sinusoidally  loading  between 

44.5  and  534  N  with  a  mean  load  of  245  N.  The  total  number  of 

cycles  was  10  .  The  final  fatigue  crack  length  was  typically  0.08 

cm.  Specimens  were  fatigued  after  hydrogen  charging  since 

niobium-hydrogen  alloys  fatigue  more  rapidly  using  low  frequency 

(35) 


loading  rates  than  at  high  frequency  loading  rates. 


I 
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TABLE  4 

Typical  Interstitial  Concentrations  for  Charged 
Niobium  Double  Torsion  Specimens 
(Concentrations  in  at.  %) 


Specimen 

Annealed 

H 

0 

N 

C 

14 

Yes 

0.55 

0.035 

0.016 

0.015 

19 

No 

2.3 

0.056 

0.031 

0.006 

49 


No 


0.30  0.068 


0.018 


0.009 
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To  determine  any  discrepancy  between  the  theoretical  and  measured 
values  of  for  niobium  d.t.  specimens,  the  compliance  as  a  function 
of  crack  length  was  measured.  Incremental  crack  lengths  were  cut  with 
a  jeweler's  saw  at  an  inclination  which  approximated  the  actual  crack 
front  profile.  Results  of  the  compliance  measurements  at  77  K  and 
273  K  are  presented  and  discussed  in  section  5.2.1. 


Single  edge  notch  specimens  were  tested  in  a  tensile  cold 

stage  designed  to  be  used  in  a  JEOL  JSM-U3  scanning  electron 

(2) 

microscope.  Load  was  controlled  by  a  variable  speed  motor  and 
measured  by  a  piezoelectric  quartz  load  cell.  Strain  and  strain 
rate  were  recorded  by  gear  driven  potentiometers.  A  video  tape 
recorder  monitored  the  propagating  crack  by  recording  the 


video  signal  from  the  s.e.m. 


The  basic  equation  that  was  used  for  calculating  for 
the  s.e.n.  specimen  has  been  given  by  Gerberich  and  Katzv  as 

K  »  (P/hW)a1/2f(a/W)  (22) 

where  P  *  applied  load, 

h  »  specimen  thickness, 

W  *  specimen  width, 

a  «  crack  length  measured  from  specimen  edge, 
f(a/'W)  =  2 (sin“(Ta/2W)  +  sec2 ( Ta / 2W) ] .  (23) 

This  expression  agrees  within  1’;  of  Brown  and  Srawley 1  s  '  fourth 

order  polynomial  for  for  values  of  ( a / V )  less  than  0.65. 

For  s.e.n.  experiments  used  in  this  research,  the  notch 

was  ”u"  shaped  to  serve  as  a  stress  concentrator  3nd  Eqn.  (22)  was 

slightly  modified  to  take  this  into  account.  The  modification  was 

in  the  form  of  a  multiplicative  factor  obtained  bv  the  graphical 

(  38 ) 

solutions  of  Peterson  for  the  effect  of  stress  concentrators 

s.e.n.  specimens.  As  an  illustration  for  the  tvpical  amount 
correction  required  to  Eqn.  (22),  f(a/W)  versus  (a/W)  is  shown  in 
Fig.  6  for  s.e.n.  specimen  119.  The  dashed  portion  of  the  curve 
represents  the  effect  of  the  "u"  shaped  notch  on  f(a/W).  Data 
acquisition  was  normallv  restricted  to  (a/W)  greater  than  0.15  and 
less  than  0.70.  For  (a/W)  greater  than  0.70,  f(a/W)  becomes 
increasingly  less  accurate.  Difficulties  wiht  the  s.e.n.  specimen 
tests  will  be  discussed  in  the  Results  Section. 
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4.3.2  Doable  Torsion  Specimens 

4 . 3 . 2 . 1  Constant  Deflection  Tests  (Decreasing  K^.) 

An  Instron  loading  machine,  floor  model  F/CM,  was  used  for 
all  the  d.t.  tests.  Figure  7  shows  the  d.t.  loading  stage  which 
was  attached  to  the  bottom  of  the  Instron' s  moveable  crosshead.  A 
typical  loading  procedure  for  constant  deflection  would  consist  of 
the  following.  A  charged  specimen  was  rapidly  loaded  to  a 
preselected  value  (loading  rate,  r^,  was  controlled  by  the  crosshead 
speed  which  was  0.2  ctn/min) .  Then  the  crosshead  was  stopped  and 
the  Instron's  load  cell  and  a  ^inear  variable  differential  transformer 
(which  measured  crosshead  displacement)  were  interfaced  with  a  Wang 
720C  computer  that  calculated  and  printed  out  the  crack,  velocity, 
stress  intensity  factor  and  the  error  associated  with  these  parameters. 
Specimens  temperatures  were  controlled  by  use  of  inert  freezing  mixture 
baths  or  by  a  regulated  flow  of  nitrogen  gas. 

4. 3.2.2  Plastic  Correction  for  Constant 
Deflection  Tests 

The  equations  given  in  the  Theory  section  apply  strictly  for  the 
case  of  a  crack  propagating  in  linear  elastic  medium.  For  niobium  this 
is  rarely  the  situation  except  at  verv  low  temperatures  near  '7  K. 
Therefore,  at  temperatures  above  77  K,  there  exists  two  contributions 
to  the  load  relaxation  in  a  hvdrogen  charged  specimen  being  fractured 
under  load.  One  contribution  is  from  the  propagating  crack.  The 
other  is  due  to  plastic  deformation  caused  by  low  temperature  creep. 
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In  order  to  correct  for  this  plasticity  contribution,  the  procedure 
outlined  below  was  followed. 

1.  A  load  relaxation  curve  using  a  hydrogen  charged 
specimen  was  recorded  on  a  strip  chart,  as 
described  in  section  4. 3.2.1. 

2.  The  test  was  repeated  using  an  uncharged  specimen  under 
the  identical  conditions  and  similarly  recorded. 

3.  Strip  charts  from  steps  1  and  2  were  superimposed  as 
shown  schematically  in  Fig.  8. 

4.  At  each  point  in  time,  t,  the  load  relaxation,  P  (t), 
due  only  to  the  advancing  crack,  tip  was  given  by 

P  (t)  =  P  (t)  -  P  (t)  +  P  (24) 

c  m  p  o 

where  P  (t)  =  load  relaxation  curve  from  the  charged  specimen, 
m 

P  (t)  ■  load  relaxation  curve  from  the  uncharged  specimen, 

P 

P  *  initial  load  both  the  uncharged  and  charged 
0  specimens  are  brought  to  at  the  moment  the  crosshead 
is  stopped  (t  =  0) . 

5.  Time  versus  P  (t)  was  plotted  and  using  Eqns.  (14')  and 
(7’)  the  corrected  velocities  and  stress  intensity 
factors,  respectively  were  calculated. 

6.  Velocity  versus  stress  intensity  data  points  from  step  5 
were  plotted. 

Shown  in  Fig.  8  is  a  region  of  the  superimposed  curves  which 
are  in  complete  coincidence  and  is  labeled  "t  ".  This  appears  to 
represent  a  region  where  no  c racking  occurs  in  the  charged  specimen 
and  can  be  explained  by  assuming  the  existence  of  an  incubation  time 
for  a  growing  hydride  to  reach  a  critical  size  before  cracking.  A 
systematic  study  of  this  incubation  period  was  not  undertaken.  However, 
it  was  found  that  the  incubation  period  was  longest  at  the  low 
temperatures  and  shortest  at  high  temperatures  with  times  ranging 
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from  a  few  seconds  to  over  a  minute.  A  theoretical  estimate  of  the 
critical  hydride  size  is  given  in  the  Results  sections. 

4. 3. 2. 3  Constant  Load  Tests  (Constant  KT) 

—  x 

The  procedure  for  constant  load  testing  was  analogous  to  the 
constant  deflection  method.  However;  in  this  case  a  constant  load  was 
maintained  instead  of  constant  deflection.  This  was  accomplished 
by  means  of  an  Instron  servo  load  control  system  which  will 
automatically  maintain  a  desired  load  by  changing  the  crosshead 
position. 

Change  in  crosshead  displacement  was  monitored  as  a  function 
of  time  and  the  calculated  crack  velocities  along  with  its  errors  were 
printed  out  on  a  teletype.  It  w^s  instructive  to  note  that  the 
constant  load  tests  do  not  yield  as  much  data  from  each  d.t.  specimen 
compared  to  the  constant  deflection  test.  One  d.t.  specimen  usually 
generates  an  entire  (V-K)  curve  using  the  constant  deflection  procedure 
whereas  only  one  or  two  data  points  are  obtainable  using  the  constant 
load  method.  Many  specimens,  therefore,  may  De  required  to  complete  a 
(V-K)  curve  for  each  temperature  and  hydrogen  concentration  when  using 
the  constant  load  method. 

4 . 3.2.4  Plastic  Correction  for  Constant  Load  Tests 

Two  contributions  to  the  change  in  crosshead  displacement  as 
a  function  of  time  will  occur  and  these  are  analogous  to  the  contri¬ 
butions  to  the  load  change  in  the  constant  deflection  tests.  The 
constant  load  test  will  result  in  greater  plasticity  corrections 
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than  the  constant  deflection  test  for  equal  hydrogen  concentrations, 
test  temperatures  and  stress  intensities. 

Fig.  9  schematically  shows  the  crosshead  displacement  as  a 
function  of  time  for  charged  and  uncharged  specimens  as  recorded  on 
a  strip  chart.  The  correction  procedure  used  has  been  summariEed 
below . 

1.  Crosshead  displacement  was  recorded  on  a  strip  chart 
for  a  hydrogen  charged  specimen. 

2.  Step  1  was  repeated  using  an  uncharged  specimen. 

3.  Using  Eqn.  (18'),  Che  crack,  velocities  V  for  the 

charged  specimen,  and  V  for  the  uncharged  specimen 
were  calculated.  p 

4.  The  differences  between  V  and  V  at  equal  times 
represents  the  velocity  dBe  onlymto  the  advancing 
crack,  V c,  i.e. 

Vc  =  VC)  -  Vp(t)  •  (23) 

The  velocity  due  only  to  the  advancing  crack  tip  was  found  to 
be  essentially  a  constant,  independent  of  time  and  represented  one 
data  point  on  a  (V-K)  curve.  An  additional  test  was  required  for 
every  change  in  load,  i.e.  stress  intensity  and/or  temperature. 


5.  RESULTS 


The  accuracy  of  measurements  for  crack  velocities  and 
stress  intensities  using  s.e.n.  and  d.t.  specimens  were  limited 
by  various  uncertainties  in  the  parameters  of  Eqns.  (7),  (14), 

(18),  (19)  and  (22).  By  estimating  these  uncertainties,  conservative 
probable  errors  were  calculated  for  V  and  K^.  The  error  analyses 
are  presented  in  Appendix  C.  Factors  limiting  the  measurements 
for  maximum  and  minimum  crack  velocities  using  the  s.e.n.  and  d.t. 
specimens  are  discussed  in  Appendix  B. 

5 . 1  Single  Edge  Motch  Data  (Increasing  K.^) 

Use  of  the  scanning  electron  microscope  tensile  stage,  though 
initially  successful,  met  with  experimental  difficulties  on 
attempting  to  examine  a  wide  range  of  crack  velocities,  hydrogen 
concentrations  and  temperatures.  Desired  strain  rates  were  difficult 
to  achieve  and/or  maintain  due  to  inherent  limitations  of  the  tensile 
stage  motor.  Temperature  range  was  limited  by  the  very  slow  cooling 
rate  of  the  stage  and  thermal  losses  to  the  microscope  column. 

The  major  obstacle  encountered  using  s.e.n.  specimens  was  the 
occurrence  of  fracture  by  initiation  of  multiple  cracks.  Multiple 
crack  fronts  made  the  calculations  of  stress  intensity  factors 


difficult . 
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5.1.1  Increasing  (V-K)  Curves 

Crack  propagation  data  from  the  s.e.n.  specimens  represent 
increasing  K^  tests.  As  the  crack  advanced,  K^  increased  according 
to  the  expression  given  by  Eqn.  (21).  Table  5  lists  the  conditions 
under  whcih  the  s.e.n.  specimens  were  tested  and  includes  the  hydrogen 
concentration  in  solid  solution  at  the  test  temperature,  hydrogen 
concentration  in  hydride  form  prior  to  loading,  true  strain  to 
failure,  e^,  average  strain  rate,  <e>=e^/At  ,  stress  free  solvus 
and  test  temperature.  Figures  10-13  shows  (V-K)  results  for  these 
tests.  Curves  were  not  fitted  for  Figs.  12  and  13  due  to  a  lack  of 
data  points. 

All  s.e.n.  specimens,  except  specimen  114,  were  tested 

slightly  below  their  solvus  temperatures.  Specimens  119  and  121, 

which  were  tested  under  very  similar  conditions,  exhibited  similar 

(V-K)  behavior.  The  stress  intensity  factor  threshold  occurred 
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between  7  and  12  MPa*m  for  specimens  119  and  121,  respectively. 

Stage  II  behavior  is  evident  in  the  (V-K)  curve  of  specimen  119 
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shown  in  Fig.  10  for  values  greater  than  30  MPa*m  .  The  (V-K) 

curve  for  specimen  121,  shown  in  Fig.  11,  displays  a  similar  trend 
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in  V  as  K^  increases  beyond  30  MPavn  .  The  (V-K)  data  for  specimen 

112  shown  in  Fig.  13  suggests  the  existence  of  a  stress  intensity 
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factor  threshold  at  approximately  23  MPa*m  and  a  comparison  with 
the  data  for  specimens  119  and  121  suggests  and  increase  of  Kj.q  as 
the  test  temperature  decreases.  While  the  curves  are  not  completely 
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Velocity  (m/s) 


Stress  Intensity  (MPa^/m") 


Fig.  10.  Crack  velocity  versus  K.  for  s.e.n.  specimen  115 
200K  and  H/Nb  *  0.5%.  1 
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11.  Crack  velocity  versus  K.  for  s.e.n.  specimen  121  tested 
190  K  and  H/Nb  -  0.5%. 


Crack  Velocity  (m 


Fi*.  12.  Crack  velocity  versus  K..  for  s.e.n.  specimen  114  tested  at 
183  K  and  H/Nb  -  0.178%. 


defined  bv  Che  available  data,  chev  do  exhibit  Stage  1  and  11  behav; 
with  V  in  Scage  II  decreasing  as  the  test  temperature  decreases. 


A  detailed  examination  of  the  fracture  surfaces  for  nvdr  -gen 
embrittled  s.e.n.  specimens  has  been  presented  bv  '.r  ossbeck . 

Fracture  surface  features  for  s.e.n.  tests  in  this  studv  are  shown 
in  Figs.  14-19  and  are  in  excellent  agreement  with  his  result-..  .ill 
f  the  fracture  surfaces  ac  low  strain  rates  exhibited  brittle 
transgranular  cleavage.  Ac  higher  strain  rates,  ductile  necking 
occurred  as  the  crack  tip  outran  the  hydrogen  diffusing  t  'ward  it. 
Fig.  14  shows  the  fracture  surface  for  specimen  11^  fcr  wh:  n  fie 
crack  propagated  from  left  to  right.  Fig.  15  is  a  higher  magnificat 
'f  the  right  side  of  Fig.  14.  The  existence  of  splinters  :>n  fa- 
fracture  surface  are  quite  prominent  and  thev  have  been  positive!" 
identified  as  niob ium-nvdride  bv  electron  diffraction  unaivsi,. 

Fig.  lb  is  the  fracture  surface  if  specimen  !  ]  u  showing  the 
reference  lines  on  the  specimen’s  surface  as  well  as  snlinteri  r  •  .< 
fracture  surface.  The  presence  -f  grain  boundaries  on  trie--  f  ,  ■  .r- 
surfaces  is  verv  clear  as  shown  ln  Fig.  1 '  which  is  a  :ugher  magr ; • . 
tion  of  the  grain  boundarv  appear  ing  in  Fig.  if  Simf  ir  :•  ,  •  :r..-  - 
are  shown  in  Fig.  IS  which  shows  part  of  the  Tract  ire  *  :  i  •-  : 
specimen  121  and  Fig.  19  whi'h  sh  sgs  the  ms.ribod  refer*--,  e  . 
the  surface  of  specimen  121  after  fracture.  It  is  leat  ; • a: 
lines  aad  no  eifect  on  fie  rraragatmg  ra  k. 
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Fig.  16.  Fracture  surface  of  s.e.n.  specimen  114.  Inscribed  reference 
lines  are  visible  on  specimen's  top  surface.  Crack  propagated 
from  right  to  left. 


Fig.  17.  Higher  magnification  of  grain  boundary  appearing  in  center  of 
Fig.  16  showing  long  horizontal  splinter.  River  lines  occur 
everywhere  on  the  fracture  surface. 
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5 . 2  Double  Torsion  Data 

5.2.1  Compliance  Measurements 

Figures  20  and  21  show  the  results  of  the  conroliance  measurements 
at  77  and  273  K,  respectively.  The  linearity  of  the  measured 
compliance  data  is  further  evidence  for  the  validity  of  the  analytical 
relationships  derived  in  the  Theory  section.  At  both  temperatures, 
however,  it  can  be  seen  that  the  theoretical  values  for  (dC/da)  are 
greater  than  the  experimentally  determinated  values  and  do  not 
extrapolate  to  the  origin.  Using  Eqn.  (17),  the  ratios  of  the 
theoretical  and  experimental  values  can  be  calculated.  At  77  K,  it 
was  found  chat 


(K  ,  /K  )_  77  =  1.14 

theory  meas.  T  =  77 


and  at  273  K, 


(26) 


(K  ,  /K.  )_  ,,,  =  1.09 

theory  meas.  T  =  273 


(27) 


Ideally,  these  ratios  should  be  unity.  To  correct  for  these  slight 

discreoancies  between  K  ,  and  K  ,  an  average  correction  factor  of 

cheory  meas. 

(1.14  +  1.09) /2  =  1.12  (28) 


was  used  to  multiply  the  measured  K  values  for  all  d.t.  tests. 

The  least  squares  line  for  the  measured  compliance  was  seen  to 

( ">  3 ) 

have  a  positive  y-intercept.  This  has  been  observed  in  steel  “  and 

( 39 ) 

polymethylmethacrylate  and  was  found  to  be  a  function  of  the 


Fig.  20.  Double  torsion  specimen  compliance  as  a  function  of  crack 
length  for  calculated  and  measured  values  at  77  K.  Solid 
line  is  least  squares  fit  for  measured  data. 
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Fig.  21.  Double  torsion  specimen  compliance  as  a  function  of  crack 
length  for  calculated  and  measured  values  at  273  K.  Solid 
line  is  least  squares  fit  for  measured  data. 
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specimen's  elastic  properties  and  dimensions  but  independent  of  the 
crack  length.  This  is  of  little  consequence  so  long  as  the 

initial  notch  length  is  greater  than  half  the  specimen's  width. 

Crack  lengths  are  measured  from  the  edge  and  since  an  initial 
notch  length  of  3  cm  was  machined  into  each  d.t.  specimen,  only  data 
for  crack  lengths  greater  than  this  minimum  value  are  shown. 

5.2.2  Constant  Deflection  (Decreasing  K^-) 

The  constant  deflection  method  represents  decreasing  tests. 
As  derived  in  the  Theory  section,  the  initial  load  is  proportional 
to  the  maximum  K  for  each  test  and  subsequent  load  relaxations 
corresponds  to  the  decreasing  Kj.  These  results  are  divided  into 
high,  low  and  intermediate  temperature  ranges  which  must  be  related 
to  the  specimen  solvus.  The  high  temperature  range  for  the  specimen's 
studied  is  between  220  K  and  303  K.  Specimens  in  this  range  have 
been  tested  at  or  above  their  solvus.  The  low  temperature  range 
is  between  77  K  and  148  K  which  corresponds  to  specimens  tested  below 
their  solvus  and  having  a  very  small  amount  of  hydrogen  remaining  in 
solid  solution.  Specimens  tested  in  the  intermediate  temperature 
range,  148  K  to  220  K,  were  slightly  above  or  below  their  solvus  but 
still  had  a  significant  amount  of  hydrogen  in  solid  solution.  Table 
6  lists  the  conditions  for  these  decreasing  tests. 

5. 2. 2.1  High  Temperature  fV-K)  Curves 

Figure  22  shows  the  results  for  several  tests  in  the  nigh 

.  values  achieved  during  the  te‘  .s. 


i 


temperature  region.  Maximum  K 


TABLE  6 
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Fig.  22  Crack  velocity  versus  curves  for  d.t.,  high 

temperature,  decreasing  tests. 


Velocity  (m/s) 


r,  f  j 
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m  .  These  tests  were  characterizes 

and  while  the  data  suggests  a  thr<.<.- 

II  extended  only  over  a  verv  limited 

stage  character  of  the  curves  are  m-'St 

273  K.  which  were  taken  slightlv  above 

their  respective  solvus  temperatures.  The  threshold  values 

decreased  while  the  Stage  II  velocities  increased  3s  the  temperature 

increased.  The  243  K  data  showed  a  further  increase  in  K.  but  these 

lo 

data  taken  at  the  solvus  temperature  did  not  exhibit  a  three  stage 
curve.  The  stress  intensity  threshold  varied  from  8  to  18  MPa  ■  m^ 
over  the  temperature  range  303  to  243  K  and  the  maximum  veloeitv 
measured  was  4  x  10  14  m/s. 


ranged  from  22.5  to  27.5  MPa  • 
by  a  strong  dependence  of  V  on 
stage  curve  the  constant  Scage 
range  of  values.  The  three 
evident  in  the  data  at  303  and 


5. 2. 2. 2  High  Temperature  Ftactographv 

Fracture  surfaces  for  the  high  temperature  tests  are  character¬ 
ized  by  predominant  brittle  failure  accompanied  by  some  ductile 
tearing  as  shown  in  Figs.  23  and  24  where  significant  amounts 
of  secondary  cracking  can  also  be  observed.  The  occurrence  of  ductile 
tearing  did  not  vary  appreciably  with  K^.  These  features  occur  along 
the  entire  fracture  surfaces. 


5. 2. 2. 3  Low  Temperature  (V-K  Curves 


Tests  at  77  K  for  the  low  temperature  region  are  shown  in 

25.  At  77  K,  the  data  are  characterized  bv  crack  velocities  whicn 

even  more  strongiv  dependent  on  K.  than  the  high  temperature  i  V— K 

i 

data.  The  values  attainable  are  a  factor  of  two  to  three  times 


ire 


greater  chan  those  tor  the  high  temperatures.  Data  taken  at  K  •» : :  •. 

the  decreasing  K^.  tests  appear  to  depend  on  the  initial  K,  value- 

for  each  test,  i.e.  the  (V-K)  data  was  multivalued.  As  seen  in 

Fig.  25,  the  velocities  observed  lay  in  the  range  10  3  to  10  ms 
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tor  Kj  values  between  30  and  70  MPa  •  m  '  and  tormed  three  separate 
curves  for  the  three  tests  carried  out.  In  this  temperature  range  tt 
is  evident  that  (V-K)  data  may  be  compared  only  when  tests  (and 
particularly  the  initial  KT  values)  are  carried  out  under  identical 
conditions.  All  three  77  K  (V-K)  curves  had  identical  slopes.  The 
dependence  of  the  data  on  the  initial  K^  value  precludes  any  discussion 
of  threshold  stress  intensities. 

Figure  26  shows  results  comparing  (V-K)  curves  for  deuterium  ir.d 
hydrogen  charged  specimen  tested  under  identical  conditions  and  an 
initial  K^  value  of  32.5  MPa  •  m^^  at  77  K.  The  curves  are  functional! 
identical  except  the  crack  velocities  from  the  deuterium  charged 
specimen  are  a  factor  of  approximatelv  15  less  than  the  crack 
velocities  of  the  hvdrogen  charged  specimen  at  equal  values. 
Comparison  with  the  expected  crack  velocity  differences  between 
hvdrogen  and  deuterium  charged  specimens  will  be  presented  in  the 
Discussion  section. 

The  results  shown  for  112  K  in  Fig.  25  are  significantly 
different  than  for  77  K.  Stress  intensity  values  less  than  Co  MPa  •  m‘ 
for  the  112  K  data  are  unavailable  due  to  the  specimen's  failure  at 
that  point.  The  112  K  (V-K)  curave  resembles  the  (V-K)  behavior  of 
the  intermediate  temperature  tests  to  be  presented  shortly. 


77  K,  Ts  -  193  K 
I  iZ  K,  Ts  =  208  K 
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Low  Temperature  Frac tographv 
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All  low  temperature  specimens  exhibited  crystallographic 
transgranular  cleavage  for  all  values  of  K^.  Figures  27  and  28  show 
the  fracture  surface  for  the  deuterium  charged  specimen.  Secondary 
cracking  was  generally  not  observed.  Presence  of  some  microvoid 
coalescence  as  well  as  river  patterns  can  be  observed  as  is  consistent 
with  the  low  hydrogen  diffusivity  at  77  K.  Hydrogen  charged 
specimens  tested  at  77  K  exhibit  identical  fracture  surface  features 
as  the  deuterium  charged  specimens  as  shown  in  Figs.  29  and  30. 

5. 2. 2. 5  Intermediate  Temperature  (V-K)  Curves 

Crack  velocity  versus  Kj.  curves  for  intermediate  temperatures 
shown  in  Figs.  31-35  are  strikingly  different  from  both  the  high  and 
low  temperature  data.  In  general,  the  (V-K)  dependence  is  not 
monotonic  but  exhibits  a  decrease  in  V  and  K^.  decreases  followed  by  an 
increase  to  a  maximum  and  a  further  decrease  in  V  with  decreasing  K . 
values.  This  behavior  is  observed  both  above  and  below  the  solvus 
temperatures  in  this  intermediate  temperature  range. 

Figure  31  shows  the  (V-K)  curves  for  218  K  at  H/Nb  =  0.46  and 
^coustic  ^mission  intensity  for  the  propagating  crack  versus  . 
Acoustic  emission  intensity  is  seen  to  be  almost  directly  proportional 
to  V.  Acoustic  emission  from  an  uncharged  specimen  used  in  the  plas¬ 
ticity  correction  •'section  4. 3. 2. 2)  was  no  greater  than  300  counts 'min 
which  is  several  orders  of  magnitude  less  to, an  a.e.  intensities  from 
the  propagating  crack.  Fig.  32  is  a  simplified  presentation  of  the 
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Fig.  27.  Fracture'  .surface  of  deuterium  oh.-  rged  d.  t.  spec  I  iron  fractured 
at  77  K.  Pronounced  crystallographic  cleavage  is  visible. 
Area  below  the  arrows  is  the  surface  of  the  machined  groove 
on  the  bottom  of  the  specimen. 
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velocity  versus  K,  of  Fig.  dl  shoving,  throe  region  behavior  of 
sing  K  tests  for  intermediate  temperatures. 


(V-K)  curve  in  Fig.  31.  Three  regions  are  labeled  A,  B  and  C  which 
are  characteristic  of  the  intermediate  temperature  results.  Since 
these  results  are  from  decreasing  K^.  tests,  Region  A  is  observed 
immediately  after  the  Instron  crosshead  is  stopped  and  is  charac¬ 
terized  by  a  decrease  in  V  to  a  local  minimum.  Specimen  53  for  which 
the  crack  was  initiated  at  the  highest  stress  intensity,  Fig.  32, 
shows  an  additional  feature  of  an  initial  increase  in  V  before  the 
local  minimum  occurs.  This  initial  increase  in  V  should  occur  for 
all  decreasing  tests  since  the  initial  V  must  start  at  zero.  With 
the  exception  of  specimen  53,  the  plastic  correction  procedure 
described  in  section  4. 3. 2. 2  prevented  the  initial  V  to  be  resolved 
in  the  other  decreasing  tests.  Region  B  corresponds  to  an  increase 
in  V  with  decreasing  K^..  Crack  velocity  increase  was  measured  to  range 
from  a  factor  of  three  to  almost  two  orders  of  magnitude.  Region  C 
is  characterized  by  a  sharp  decrease  in  V  with  decreasing  where 
V  continues  to  decrease  to  values  beyond  the  limit  of  measurement. 

Figure  33  shows  the  (V-K)  curves  for  two  specimens  tested  it  i: 
below  their  solvus  temperatures  and  having  identical  hydrogen  concen¬ 
trations.  Reproducibility  of  data  at  low  values  is  shown  to  be 
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excellent.  Crack  initiation  occurred  at  30  and  40  MPa  •  m  for  tnese 
two  specimens.  During  Region  B,  V  increased  by  almost  two  orders  ■>( 
magnitude.  Figure  34  shows  the  CV-K)  curves  for  two  specimens  having 
identical  hydrogen  concentrations  and  tested  at  either  13  or  13  K 
above  their  ,olvus  temperatures.  Regions  A  and  C  in  this  figure  are 
^uite  similar  to  Regions  A  and  C  of  Figs.  32  and  33.  Region  B  of 
Fig.  34,  however,  shows  onlv  an  increase  of  a  factor  of  three  in  V 


before  Region  C  is  reached.  Figure  35  shows  all  the  intermedi ite 
temperature  data  plotted  on  one  graph  for  comparison. 


Fracture  surface  features  of  the  intermediate  temperature 
range  tests  reflected  the  behavior  of  the  three  region  iV-K)  curves. 
Figure  36  shows  the  notch  area  where  cracking  was  initiated  in  Keg;  m 
The  fracture  surface  was  predominantly  transgranular  cleavage  with  a 
substantial  amount  of  secondary  cracking.  The  area  on  the  fracture 
surface  near  the  top  surface  of  the  specimen  exhibits  some  ductile 
shearing.  Figure  37  is  the  center  area  of  Fig.  36  at  higher  m.rng  i :  -.a 
tion  and  shows  the  presence  of  brittle  fracture  along  with  some 
ductility.  The  bottom  area  of  the  fracture  surface,  near  the  groove, 
shows  only  brittle  failure  with  a  multitude  of  river  lines  as  shown 
in  Fig.  38. 

Fracture  surfaces  corresponding  to  Region  3  behavior  ••xhihited 
predominantly  cleavage  with  a  greater  amount  of  ductile  fracture  ‘-nun 
Region  A  as  shown  in  Fig.  39.  The  center  and  bottom  area  of  the 
fracture  surface  corses  lag  to  Region  3  also  shows  i  sign;:  :  -in', 
amount  of  seconds.  ciacki.  „.  Figures  ->0  ind  -1  show  t  t  ran  s . :  :  n 
region  between  ductile  and  Hrittle  fracture  .n  Region  B  it  hie:,  ma i- 
nificacion.  These  figures  also  show  the  existence  of  small  s',  ip 
steps  m  the  rumpled  shear  surfaces.  The  ductile  shear  ire  i  v  is  n  t 
entire  Lv  without  brittle  cracks.  Figures  -3  o.u  •  >  s-..-v  mo-  . 

•'  arres ponding  to  Region  3  wnich  clear.  .-  shows  tne  presence  o:  '  . 
dir  tile  mu'  brittle  fracture  i  long  wira  sc  ruarv  rue  r'. :  n  g . 


Fig.  38.  Fracture  surface  resulting  from  Region  A  behavior  near  the 
bottom  groove  (below  arrows).  Predominance  of  river  lines 
are  visible. 


Fig.  39.  Fracture  surface  resulting  from  Region  B  behavior  showing 
top  of  d.t.  specimen.  Significant  amount  of  ductile 
tearing  can  be  seen  as  well  as  secondarv  cracking. 
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Region  C  behavior  was  typically  characterized  by  brittle 
cleavage  and  the  presence  of  river  as  shown  in  Fig.  44.  In  contrast 
to  Regions  A  and  B  little  evidence  for  ductile  fracture  was  observed. 
Also  visible  in  Fig.  44  are  elongated  grain  boundaries  along  the 
rolling  direction  of  the  plate. 

5.2.3  Constant  Load  (Constant  K^.)  Data 

These  results  represent  (V-K)  data  obtained  from  constant 
stress  intensity  tests  as  described  in  section  4. 3. 2. 3.  Crack 
propagation  for  constant  tests  was  significantly  different  from  the 
decreasing  tests.  The  constant  data  exhibited  the  familiar 
three  stage  velocities,  mentioned  in  section  2.3,  often  seen  in 
stress  corrosion  cracking. 

5.2. 3.1  Constant  KT  (V-K)  Curves 

”  '  A. 

Table  7  lists  the  d.t.  constant  test  conditions  and 
Fig.  45  shows  the  results  for  the  constant  tests  at  four  tempera¬ 
tures,  180,  218,  232  and  303  K.  At  the  three  highest  temperatures, 
tests  were  carried  out  above  the  specimen's  solvus  temperatures  while 
the  180  K  test  was  carried  out  at  42  K  below  the  specimen's  solvus. 

Onset  of  rapid  crack  growth  in  Stage  III  was  found  to  occur  at  ->]  : 

1/2 

1.2  MPa  •  m  and  independent  o  f  temperature  or  hydrogen  concent  r  it  tor. . 
n  definite  threshold  stress  intensicv  factor,  KT  ,  below  which  ni> 
crack  velocity  (within  the  limits  of  the  measurement  capab i 1  i  tv ;  .  >u!  : 
be  measured  appears  to  exist  for  all  the  temperatures  examined.  it 
value  of  Kt  decreased  as  the  test  temperature  decreased  in  ••entrust  t 


f 


i 

t 


j 


4gfi‘ 
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results  obtained  from  the  increasing  (s.e.n.  data)  and  decreasing 

1/2 

Kj.  tests.  The  threshold  1^  values  were  25  r  1  MPA  •  m  at  30)  K, 

22.3  ±  1  MPa  •  m^~  at  232  and  218  K  and  15  t  1  Mia  •  at  180  k. 

In  Stage  I,  V  was  observed  to  increase  by  several  orders  . t  magn i t ude 

1/2 

for  an  increase  in  K  by  only  4  or  5  MPa  •  m 

A  definite  Stage  II  constant  (V-K)  region  was  observed  for  ill 
four  temperacures .  The  velocities  in  Stage  II  decreased  as  the, 
temperature  decreased  from  about  10  m/s  at  303  K  to  5  x  10  m/'s  at 
180  K.  The  trend  with  temperature  is  the  same  as  that  observed  far 
increasing  K.^  (s.e.n.  data)  and  decreasing  curves  but  the  velocity 
values  are  orders  of  magnitude  lower  for  the  constant  measurements. 

5 . 2 . 3 . 2  Constant  Fractography 

As  mentioned  in  the  Experimental  Procedure  section,  all 
constant  K ^  d.t.  specimens  were  fatigued  prior  to  crack  velccitv 
measurements.  The  fatigue  surface  of  d.t.  specimen  51  and  the 
fatigue  boundary  which  is  delineated  by  a  white  line  is  shown  in 
Fig.  46.  Figures  47  and  48  shows  a  higher  magnification  of  the  f  align* 
crack  fracture  surface  which  exhibits  a  predominance  of  ductile 
tearing  as  is  expected.  Figure  49  shows  the  interface  between  the 
fatigue  crack  and  constant  KT  test  fracture  surfaces.  Brittle 
cleavage  steps  can  be  clearly  seen  on  the  side  of  the  interf.i<  *■ 
corresponding  to  the  constant  K  crack  propagation  and  ductile 
features  are  seen  on  the  fatigue  surface  or.  the  other  side. 

Fracture  surface  features  for  constant  K.  specimen.:-  test-1 
above  and  below  their  solvus  temperatures  differed  slight]-.-  fr>m  u-_ 


another.  Figure  50  shows  the  fracture  surface  near  the  notch  area  r 

d.t.  specimen  36  tested  above  its  solvus.  The  top  of  the  fracture 
surface  shows  the  presence  of  some  slight  ductile  shearing.  The 
majority  of  the  specimen's  fracture  surface  exhibited  transgranular 
cleavage  with  some  ductile  fracture,  much  secondary  cracking  and  an 
abundance  of  river  lines  due  to  brittle  cleavage  steps  which  ire  sn.-wr; 
at  higher  magnifications  in  Figs.  51  and  52.  Figures  53  and  3-  shew  the 
fracture  surface  for  the  specimen  tested  at  180  K  which  is  12  K  below 
its  solvus.  Its  fracture  surface  exhibited  similar  secondary- 
cracking  and  the  profusion  of  river  lines  similar  to  those  shown  by 
the  specimens  tested  above  their  solvus  temperatures.  The  specimen 
tested  below  its  solvus  did  not,  however,  exhibit  anv  of  the  ductile 
shearing  that  was  present  for  specimens  tested  above  their  solvus. 

5.2.1  Critical  Hydride  Sice 

As  mentioned  in  the  Experimental  Procedure  section,  an 
incubation  interval,  during  wnich  no  crack  propagation  occurred,  became 
apparent  during  the  plastic  correction  procedure  for  measuring  d.t. 
specimen  crack  velocities.  A  critical  hvdride  size,  1.,  before  crack 
advance  occurs,  can  be  estimated  bv  assuming  stress  assisted 
diffusion  of  hydrogen  into  the  crack  tip  region.  Under  these  concii:  ;  ons  . 
a  hvdride  will  continue  to  grow  until  i  critical  site  is  reached  where 
the  local  stress  distribution  exceeds  the  hvdride' s  fracture  stress 


and  failure  occurs. 
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The  assumptions  in  estimating  the  critical  hydride  size  are 
listed  below. 

1.  The  stress  induced  hydride  has  a  composition  of 
NbH  75  (Cjjyd  -  43  at  .2). 

2.  The  hydride  has  the  shape  of  an  ellipse  of  revolution 
whose  major  axis,  a,  is  twice  the  length  of  the  minor 
axis,  b. 

3.  The  surrounding  volume  from  which  hydrogen  diffuses  is 
bounded  by  a  sphere  of  radius  r  where  r  is  much  greater 
than  the  length  a. 

4.  Stress  gradient  effects  are  not  taken  into  account. 

For  a  hydride  to  form  in  an  alloy  with  C°  less  than  43  at  .2,  hydrogen 
must  diffuse  into  the  volume  of  the  hydride  from  the  surrounding  spher¬ 
ical  volume.  Therefore, 

0.43irab2  -  (4/3)itr3C°  (29) 

The  incubation  interval  t.  ,  for  the  critical  hydride  size  is  given 

me 

by 

r  (301 

where 

dJ*5  -  2.4  x  10~4exp(-8.20kJ/mol/RT)  cm/s  . (41)  (31) 

Substituting  Eqn.  (30)  into  Eqn.  (29)  and  letting  lc  ■  a  «  2b,  yields 

lc  -  76.7<C°)L/3  exp(-4. lkJ/mol/RT)  •  (tlnc)1/2  (32) 

where  C°  -  hydrogen  concentration  in  solution  at  test 
temperature  in  at.  2. 
lc  *  critical  hydride  length  in  um. 
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Table  8  lists  the  estimated  hydride  si2e  using  Eqn.  (32)  for  various 

test  temperatures  and  hydrogen  concentrations.  Metallographic  obser- 

(2) 

vations  of  stress  induced  hydrides  in  this  study  and  Grossbeck's 

for  niobium-hydrogen  alloys  showed  sizes  ranging  from  10  to  50 

Um,  in  agreement  with  the  estimated  sizes  for  the  higher  temperatures. 


TABLE  8 

Critical  Hydride  Size 

1 

Hydrogen  Concentration 

In  Solid  Solution  At 

Test 

Temperature 

Cinc 

(s) 

(Vim) 

Test  Temperature 
(at.%) 

(  K  ) 

0.02 

6.4  x  10"6 

77 

78 

0.6 

0.0017 

112 

30 

30 

0.46 

218 

24 

29 

0.46 

221 

21 

35 

1.3 

243 

10 

41 

2.3 

273 

6 

I 


6.  DISCUSSION 


6 . 1  Single  Edge  Notch  Data  (Increasing  Tests) 

Velocity  versus  stress  intensity  results  for  the  s.e.n.  tests 

(Figs.  10-13)  can  be  explained  assuming  stress  induced  hydrides 

precipitate  in  the  region  of  stress  concentration  produced  at  the 

(2) 

specimen's  notch.  As  the  stress  is  applied,  hydrogen  diffuses  to 

the  stress  concentration  at  the  notch  due  to  the  chemical  potential 

gradient  and  stress  induced  hydrides  form  and  cleave.  Figure  55a  is  a 

schematic  representation  of  a  (V-K)  curve  for  the  increasing  tests 

and  has  been  labeled  with  points  R,  S,  T  and  U.  Figures  55D-55e  are 

schematic  representations  of  the  conditions  occurring  at  the  crack 

tip  corresponding  to  point  R,  S,  T  and  U,  respectively.  The  dashed 

ellipses  are  the  plastic  zone  boundaries  and  the  smaller  ellipses  within 

the  boundaries  are  the  precipitated  hydrides.  These  figures  will 

be  referenced  in  the  following  discussion. 

At  low  stresses  hydrogen  flux  to  the  crack  tip  takes  place  in 

response  to  the  chemical  potential  gradient  caused  by  the  external 

stress.  At  temperatures  below  the  solvus  temperature  the  stress  will 

(9  11  42) 

also  induce  hydride  precipitation  at  the  crack  tip  ’  ’  while  at 

temperatures  above  the  solvus  temperature  this  will  occur  only  at 

stresses  large  enough  to  increase  the  local  solvus  temperature  to  a 

(9  11  42) 

value  above  the  test  temperature.  ’  ’  At  a  particular  stress, 

corresponding  to  a  threshold  stress  intensity,  Kjq,  the  hydrides  will 
fracture  and  the  crack  will  propagate  to  the  interface  with  the  solid 

* 

\ 

< 


Schematic  representation  of  the  (V-K)  curve  for  increasing 
Kj  tests  (a),  and  conditions  of  tne  plastic  zone  and 
hydrides  at  the  crack  tip  as  crack  growth  occurs  (b) ,  (c), 
(d)  and  (e). 
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solution  (Fig.  55b).  The  crack  front  of  the  fractured  hydride  will 


then  serve  as  a  stress  concentrator  and  the  process  will  repeat 
itself.  (9*U'42> 


As  shown  previously  by  Kqn.  (7 2),  Kj  for  the  s.e.n.  specimens, 

increases  as  the  crack  extends.  The  increasing  would  result  in  a 

larger  plastic  zone,  a  higher  local  concentration  of  hydrogen  and  a 

decrease  in  the  local  solvus  thus  resulting  in  faster  formation  of 

hydrides  (Fig.  55c),  and  consequently  an  increase  in  crack  velocity. 

In  the  low  K  ,  increasing  velocity  region  the  stress  induced  hydride 

extends  throughout  the  plastic  zone  and  the  crack  velocity  will 

increase  as  and  the  plastic  zone  size  increases.  At  sufficiently 

high  Kj  (Fig.  55d) ,  the  growth  of  the  hydride  and  hence  the  crack 

velocity  will  become  rate  limited  by  diffusion  of  hydrogen  from  the 

surrounding  area  to  the  plastic  zone.  With  a  further  increase  in  K^, 

the  plastic  zone  size,  r  ,  will  continue  to  grow  (Fig.  55e) ,  since  r^ 

2 

is  proportional  to  K^,  but  the  crack  velocity  will  not  increase  since 

it  is  now  rate  limited  by  diffusion  of  hydrogen.  In  this  region  the 

extent  of  the  hydride  field  is  less  than  r^.  This  will  result  in  Stage 

II  behavior  of  the  (V-K)  curve.  This  is,  indeed,  what  the  s.e.n. 

(V-K)  results  showed  (Figs.  10-13). 

The  fractographic  observations  of  this  study  and  those 

(2) 

previously  reported  support  this  model.  In  particular,  a  stress 
induced  hydride  embrittlement  mechanism  can  account  for  the  observations 
that  all  fracture  surfaces  of  the  s.e.n.  specimens  exhibited  brittle 
transgranular  cleavage  at  low  strain  rates  and  ductile  necking  at  high 


strain  rates.  The  ductile  fracture  results  when  failure  occurs  at 
high  Kj  in  the  plastic  zone  in  front  of  the  precipitated  hydride 
zone. 


6.2  Double  Torsion  Specimen 

6.2.1  Decreasing  K^  Data 

The  (V-K)  behavior  for  the  decreasing  tests  can  be 
satisfactorily  explained  by  the  effects  of  the  stress  on  the  solvus 
and  its  effect  on  the  crack  velocity.  These  tests  were  carried  out 
above  K.^,  the  stress  intensity  for  fracture  of  the  hydrides,  and 
therefore  the  crack  velocity  reflects  the  rate  of  hydride  formation. 
Figure  56  shows  the  low  concentration  end  of  the  niobium-hydrogen 
equilibrium  phase  diagram^  and  the  relative  positions  of  the  high, 
low  and  intermediate  temperature  tests  with  respect  to  their  solvus. 

The  (V-K)  behavior  of  these  three  temperature  ranges  will  be 
discussed  below. 

Figure  57a  is  a  schematic  representation  of  the  (V-K)  curves  at 
high  temperatures.  Curves  for  two  temperatures,  T^  and  T2,  are  shown 
where  T^  is  greater  than  T2.  The  curve  for  T^  is  labeled  with  points 
X  and  Y.  Curve  T2  is  labled  with  point  Z  where  points  X  and  Z 
represent  the  same  K^  value.  Figures  57b-57d  are  schematic  representa¬ 
tions  of  the  conditions  at  the  crack  tip  for  points  X,  Y  and  Z, 
respectively.  Reference  will  be  made  to  these  figures  for  the  follow¬ 


ing  high  temperature  discussion. 


Temperature 


(  .(High  Temperature  Data) 
r  ^Solvus 


Three  Region  Behavior 
(Intermediate  Temperature  Data) 


■(Low  Temperature  Data) 
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Fig.  56.  Low  concentration  end  of  niobium-hydrogen  equilibrium  phase 
diagram  showing  positions  of  decreasing  K-  tests  relative  to 
the  solvus.  1 


(point  Y) 


(point  Z) 

/ 
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Fig.  57.  Schematic  representations  of  (V-K)  curves  for  high 

temperature  decreasing  tests  (a),  and  the  corresponding 
hydride  and  plastic  2one  conditions  at  the  crack  tip  (b) , 
(c)  and  (d) . 
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If  a  specimen  is  tested  above  its  solvus  temperature 

(T  >  Tg  +  35  K)  then  only  a  very  high  stress  should  be  able  to  stress 

induce  hydrides  at  the  crack  tip.  It  can  be  shown  that  the  local 

solvus  temperature  at  the  crack  tip  is  increased  about  25  K  for  a 

local  triaxial  stress  of  5%  of  the  shear  modulus  with  1.0  at .X  hydrogen 

(2) 

in  solution.  In  addition  to  the  effect  of  elastic  stress  fields  at 

the  crack  tip  on  raising  the  hydride  solvus  temperature,  the  plastic 
deformation  will  assist  the  formation  of  hydrides  by  reducing  the  free 
energy  associated  with  the  elastic  and  plastic  accommodation  terms 
Thus,  on  loading  the  specimen,  crack  propagation  will  not  occur  until 
the  stress  intensity  and  the  crack  tip  plasticity  raises  the  local  solvus 
temperature  to  allow  hydride  preciptation,  (Fig.  57b) .  Once  the 
hydrides  are  precipitated,  they  will  fracture  and  the  crack  will 
again  propagate  at  a  rate  which  will  be  determined  by  the  race  of 
hydride  formation  at  the  crack  tip.  Since  decreases  as  the  crack 
propagates,  the  driving  force  for  hydride  precipitation  and  corres- 
poindingly  the  rate  of  hydride  formation  will  decrease.  Figure  57c  is  a 
schematic  representation  of  the  smaller  plastic  zone  and  hydrides 
due  to  a  drop  in  resulting  in  a  decrease  in  crack  velocity.  The 
high  temperature  (V-K)  curaves,  (Fig.  22),  exhibited  this  behavior. 

For  high  temperature  tests,  the  local  solvus  temperature  on  the 
crack  tip  must  be  increased  to  the  testing  temperature,  by  the  local 
stress  field  and  plastic  deformation,  to  allow  the  hydrides  to 
precipitate  and  cleave.  This  determines  the  at  which  the  crack 
begins  to  propagate  during  loading  of  the  specimen.  The  relatively 
low  value  of  at  which  the  crack  begins  to  propagate  compared  to  the 


results  at  low  and  intermediate  temperatures,  (Figs.  25  and  35),  may 

be  explained  by  the  fact,  that  at  high  temperatures,  the  Nb  yield 

stress  is  quite  low  (^350  MPa).  This  results  in  a  relatively  low 

applied  stress  required  to  perform  the  plastic  work  which  assists  the 

plastic  accommodation  of  the  hydrides,  allowing  them  to  be  stress 

induced.  It  does  suggest  that  the  dominant  terms  in  shifting  the 

solvus  temperature  at  the  crack  tip  are  the  reduction  of  the  elastic 

(9) 

and  plastic  accommodation  free  energies.  This  is  a  conclusion  which 
is  consistent  with  previous  suggestions. 

As  can  be  seen  in  Fig.  22  the  (V-K)  relation  in  this  temperature 
range  is  not  single  valued  as  the  curves  depend  on  the  temperature  of  the 
test.  At  a  constant  K^.  value  the  velocity  obtained  at  higher 
temperature  is  greater  than  Chat  obtained  or  extrapolated  from  a  lower 
temperature.  This  is  not  unexpected  since  the  temperature  dependence 
of  the  crack  propagation  at  high  temperatures  will  reflect  the 
diffusion  of  hydrogen  to  the  crack  tip  as  well  as  any  variation  of 
crack  tip  plasticity  with  temperature  since  the  dominant  effect  of 
stress  appears  to  be  the  plastic  accommodation  terms.  The  conditions  at 
the  crack  tip  responsible  for  the  smaller  velocity  at  equal  Kj  at 
temperature  T£  are  shown  in  Fig.  57d.  The  increase  in  yield  strength 
at  the  lower  temperature  is  responsibile  for  a  smaller  plastic  zone 
and  hydride  field  resulting  in  lower  crack  velocities  as  compared  to 
the  curve  for  T^  at  equal  K^.. 

The  high  temperature  fractography  (Figs.  23  and  24)  shows  the 
presence  of  some  ductile  tearing  which  is  consistent  with  the  expected 


increase  in  plasticity  for  these  temperatures.  These  regions  of 
ductile  fracture  probably  occurred  between  the  stress  induced 
hydrides . 

Figure  58a  is  a  schematic  representation  of  the  decreasing 
Kj  (V-K)  curves  at  low  temperatures.  Curves  for  two  temperatures, 
and  T,  are  shown  where  is  greater  than  T2.  The  curves  for 
and  T2  are  labeled  with  points  M  and  N  which  represent  the  same 
value.  Figures  58b  and  58c  are  schematic  representations  of  the  condi¬ 
tions  at  the  crack  tip  for  points  M  and  N,  respectively.  Reference 
will  be  made  to  these  figures  for  the  following  low  temperature 
discussion. 

For  specimens  tested  well  below  their  solvus  temperatures 

(T  <  Tg  -  60  K) ,  most  of  the  hydrogen  in  solution  will  have  already 

been  precipitated  as  bulk  hydrides  prior  to  loading  and  application 

of  a  stress  will  result  in  only  a  very  small  amount  of  stress  induced 

hydrides  at  the  crack  tip  (Fig.  58c).  These  stress  induced  hydrides 

as  well  as  any  favorably  located  hydrides  which  precipitate  on  cooling 

will  cleave  causing  crack  propagation  until  the  crack  arrests  at  the 

hydride-solid  solution  interface.  Continued  hydride  formation  and  the 

concurrent  crack  propagation  will  proceed  slowly  since  the  hydrogen 

concentration  in  solid  solution  in  the  area  surrounding  the  crack  is 

small  and  the  diffusivity  is  low  at  low  temperatures.  Crack 

propagation  would  be  controlled  by  the  rate  of  stress  induced  hydride 

formation  as  the  tests  were  carried  out  above  KT  .  Since  these  tescs 

Io 

are  carried  out  below  the  solvus  the  rate  of  hydride  formation  is 
limited  by  the  hydrogen  diffusion  to  the  crack  tip. 


w 


Fig. 
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58.  Schematic  representations  of  (V-K)  curves  for  low  temperature 
decreasing  K  tests  (a),  and  the  corresponding  hydride 
conditions  at  the  crack  tip  (b)  and  (c) . 


In  the  low  temperature  mage  the  Nb  yield  stress  in  high 
(vl400  MPa)  and  the  crack  tip  plastic  zone  correspondingly  small 
(Fig.  58c).  In  the  absence  of  a  large  plastic  zone,  the  local  K^  will 
be  large  (a  relatively  sharp  crack),  and  the  dominant  effect  of  stress 

(9) 

on  the  hydride  formation  will  be  on  the  elastic  free  energy  terms. 

The  local  stress  field  at  the  crack  tip  provides  the  driving  force  for 
hydrogen  diffusion  and  hydride  formation.  The  rate  of  hydride 
formation  and  hence  of  crack  propagation  would  reflect  primarily  the 
hydrogen  flux  to  the  crack  tip  region  and  the  temperature  dependence  of 
the  (V-K)  curve  would  then  reflect  the  temperature  dependence  of  the 
hydrogen  diffusivitv.  For  the  same  K^  values  (Fig.  58a),  the  hydride  at 
the  crack  tip  as  well  as  the  crack  velocity  at  T^  will  be  greater  than 
the  crack  velocity  at  T^  due  to  the  lower  hydrogen  diffusivity. 

Since  K^  decreases  as  the  crack  propagates,  the  driving  force 
for  hydride  formation  at  the  crack  tip,  i.e.,  the  stress  induced  decrease 
in  the  free  energy  of  the  hydride,  correspondingly  decreases  and  the 
rate  of  hydride  formation  decreases.  This  behavior  would  correspond  to  a 
steep  decrease  in  the  crack  velocity  at  a  relatively  high  K^  and  this  is, 
in  fact,  what  is  observed  at  77  K  as  shown  previously  in  Fig.  25. 

In  this  low  temperature  range  it  is  striking  that  decreasing 
tests  begun  at  different  K^  values  at  the  same  temperature  do  not  fall 
on  a  single  (V-K)  curve.  Tests  begun  at  low  K^  values  show  much  higher 
crack  velocities  than  those  extrapolated  from  data  obtained  in  tests 
initiated  at  higher  Kj  values  (Fig.  25).  Figure  59a  is  a  schematic  repre¬ 
sentation  of  the  multivalued  nature  of  the  decreasing  K^  (V-K)  curves 
at  77  K.  The  curve  labeled  with  point  0  is  the  low  initial  K^  curve 


Fig.  59 
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(a) 


cracked  hydride 


ppt.  hydride 
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Schematic  representation  of  multivalued  (V-K)  curves  for  the 
decreasing  KT  tests  at  77  K  (a),  and  the  corresponding  h-.-dride 


conditions 
and  (c) . 


at  the  crack  tip  for  different  values  (b) 


while  Che  curve  labeled  wich  point  P  is  Che  high  ini'.ial  curve. 

The  fact  that  the  cuzves  are  not  colinear  illustrates  the  double 
valued  nature  of  the  (V-K)  behavior  at  77  K.  Figures  59b  and  59c  are 
schematic  representations  of  the  condit  ions  at  the  cra<~k  tip  for  points  0 
and  P,  respectively.  Reference  will  be  made  to  chese  figures  in  the 
following  discussion. 

The  multivalued  nature  of  the  (V-K)  data  at  77  K  can  be  under¬ 
stood  in  view  of  che  fact  that  a  significant  amount  of  the  hydrogen  was 
precipitated  as  hvdrides  on  cooling  prior  to  loading.  On  application  of 
a  stress,  at  a  stress  intensity  above  KT  ,  these  hvdrides  will  fracture 
and  the  measured  crack  veloeicv  will  be  the  result  of  crack  propagation 
through  the  precipitaced  hydrides  and  the  formation  of  stress  induced 
hvdrides  in  front  of  the  cracked  hydrides  as  well  as  In  front  of  the 

macroscopic  notch.  The  local  stress  intensity  at  the  macroscopic 
6  t  f 

crack  tip,  KT  ,  will  be  reduced  bv  this  hvdride  cracking  in  the  same 
1 

manner  as  bv  a  plastic  zone  at  the  crack  tip.  For  tests  which  begin  at 

low  K.  (Fig.  59b',  the  extent  of  this  fracture  zone  is  limited  and  crack 

propagation  occurs  bv  stress  induced  hvdride  formation  and  cleavage  *c 

£  t  f 

rates  determined  bv  hvdrogen  diffusion  to  the  crack  tip  at  K.,  K.. 

The  same  mechanism  aoplies  for  tests  begun  at  high  values  (Fig.  59c ) 
but  a  more  extensive  region  o:  hvdride  cracking  occurs  in  the  precipi¬ 
tated  nvdrides  thus  giving  an  effectively  large  crack  opening  displace- 

g  t  f 

ment  with  a  corresponding  decrease  in  t.ne  Iccai  K^'*.  The  crac* 

e  f  f 

propa^ces  at  hvdride  rornacion  rat^s  consistent  with  the  K;  )es» 

i 

chan  Kr  (applied).  At  the  hi<:iest  value  tested,  rnanv  nvdrides  will 


fracture  and  more  will  be  stress  induced  from  che  little  remaining 
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hydrogen  In  solution.  This  will  result  in  the  highest  crack  velocity 
as  was  observed.  Due  to  the  decrease  of  the  local  resulting  from 
the  large  crack  field,  the  velocity  will  be  lower  than  that  extra¬ 
polated  from  tests  at  low  initial  values.  The  identical  slopes 
of  the  77  K  (V-K)  data  which  were  obtained  at  different  initial 
values  (Fig.  25)  indicates  that  the  rate  controlling  mechanism, 

hydrogen  diffusion  to  the  crack  tip,  is  independent  of  the  initial 

ef  f 

K.j.  values  and  depends  on  the  .  Equal  crack  velocities  imply  equal 

0  £  £ 

hydrogen  fluxes  and  hence  equal  even  though  the  macroscopic  Kj 

values  may  differ  greatly. 

The  differences  between  the  crack  velocities  for  hydrogen  and 

deuterium  charged  specimens  (Fig.  26)  can  be  understood  by  considering 

the  isotope  effect  on  the  hydrogen  diffusion.  At  77  K,  the  ratio 

between  the  hydrogen  and  deuterium  diffusion  coefficients  using  data 

(43) 

from  Wipf  and  Alefield  is  given  by 

'7  *  12-8  •  03) 

H  D  T*  /  7 

(44 ) 

Extrapolation  of  data  by  Matusiewicz  and  Sirnbaum  yields  a  ratio  of 
17.  Since  the  crack  velocities  in  this  low  temperature  range  are 
determined  by  hydrogen  flux  to  the  crack  tip,  the  isotope  effect  on 
velocity  should  be  equal  to  its  effect  on  the  diffusivitv.  These  ratios 
agree  very  veil  with  the  value  of  15  measured  for  the  ratio  of  the 
velocities  between  hvdrogen  and  deuterium  charged  specimens  tested  at 
77  K.  It  is  clear  that  the  isotope  effect  accounts  for  the  observed 
differences  in  the  crack  velocities  between  hvdrogen  and  deuterium 
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charged  d.t.  specimens  and  that  the  crack  velocity  in  this  temperature 
range  is  controlled  by  hydrogen  flux  to  the  crack. 

Due  to  the  presence  of  hydrides  prior  to  loading,  a  predominance 
of  brittle  cleavage  should  occur  on  the  fracture  surface  and  this  was 
observed  to  be  the  case.  For  regions  where  the  crack  outran  the  supply 
of  hydrogen  microvoid  coalesence  should  occur  and  this  also  was 
observed  as  pointed  out  in  Section  5. 2. 2. 4. 

The  (V-K)  behavior  in  the  intermediate  temperature  range,  i.e. 
temperatures  in  the  vicinity  of  the  solvus,  is  considerably  more  complex 
than  at  high  or  low  temperatures  as  was  shown  in  Fig.  35.  This  com¬ 
plexity  seems  to  result  from  the  interplay  between  the  large  effects 
of  stress  on  the  hydride  precipitation  and  the  flux  of  hydrogen  into 
the  crack  tip  region  which  allows  this  precipitation  to  occur.  In 
this  temperature  range  the  variation  of  the  velocity  with  Kj.  was 
divided  into  three  regions  as  presented  previouslv  in  the  Results 
section  (5. 2. 2. 5)  and  shown  in  Fig.  32.  Figure  60  is  a  schematic  repre¬ 
sentation  of  the  hydrides  in  the  vicinity  of  the  crack  tip  as  the  crack 
advances  through  Regions  A,  B  and  C. 

Region  A,  which  occurs  immediately  on  loading  the  specimen  to 
initiate  the  crack,  consists  of  a  very  rapid  increase  of  the  crack 
velocity  as  decreases  with  the  crack  velocity  reaching  a  maximum 
value  and  then  decreasing.  As  seen  in  Figs.  33  and  34  the  initial 
increase  is  not  always  observed  and  seems  to  result  when  the  crack 
Is  initiated  at  high  stress  intensity.  At  lower  Kj  values  for  the 
initiation  of  the  fracture  only,  the  decreasing  part  of  the  (V-to  .urv.; 
in  Region  A  is  seen.  In  this  temperature  range,  where  the  hvdrogen 
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Fig.  60.  Schematic  representation  of  hydride  behavior  responsible  for 
three  region  (V-K)  curves  at  intermediate  temperatures  for 
decreasing  tests  as  crack  advance  occurs  (a),  (b)  and  to. 
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diffusivity  is  high,  large  hydride  precipitates  form  during  loading 
of  the  specimen  providing  a  continuous  crack  path  (Fig.  60a).  When 
these  hydrides  crack,  the  propagation  rate  will  increase  rapidly 
into  the  pre-existing  precipitated  hydride  field.  The  growth  of  these 
hydrides  during  loading  and  initial  crack  propagation,  results  in  a 
depletion  of  the  hydrogen  concentration  in  the  volume  in  front  of  the 
notch.  Thus  the  crack  velocity  quickly  decreases  as  the  c rack  grows 
out  of  the  region  of  hydrides  precipitated  during  loading.  For  cracks 
begun  at  lower  K^,  a  smaller  hydride  zone  is  formed  and  the  initial 
high  velocity  corresponds  to  the  crack  propagation  through  the  hydride 
field.  Only  the  decrease  in  velocity  due  to  the  depletion  of  the 
hydrogen  concentration  by  the  initial  precipitation  is  seen  because  at 
the  lower  values,  the  hydride  growth  is  relatively  slower. 

As  the  crack  slows  down  in  Region  A,  existing  hydrides  continue 
to  grow  and  new  ones  are  precipitated  ahead  of  the  crack  tip  as  depicted 
in  Fig.  60b.  This  occurs  due  to  stress  assisted  diffusion  of 
hydrogen  to  the  plastic  zone  area.  The  further  growth  of  the  hydrides 
allows  the  crack  to  accelerate  and  the  crack  propagates  by  linking  up 
of  the  cracked  hydrides  as  in  Region  B. 

In  Region  B,  the  velocity  is  determined  by  the  rate  of  hydride 
growth  and  is  controlled  by  diffusion  of  hydrogen  in  the  stress 
field  at  the  crack  tip.  The  velocity  increase  results  because  the 
hydrides  precipitate  well  in  front  of  the  crack  tip  and  the  crack  can 
advance  into  an  existing  hydride  field  (Fig.  60b).  Due  to  the  high 
crack  velocity,  the  linkage  of  the  cracked  hydrides  by  additional 
hydrides  formed  during  the  cracking  process,  may  not  be  complete  and 
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some  ductile  tearing  may  result  between  the  cleaved  hydrides  as  is 
observed  (Figs.  42  and  43). 

As  a  result  of  the  growth  of  the  hydrides  in  Regions  A  and  B, 
che  volume  in  front  of  the  crack  becomes  depleted  of  hydrogen 
requiring  diffusion  from  greater  distances  to  form  hydrides.  This 
results  in  a  decrease  of  the  hydrogen  flux.  In  addition,  the  driving 
force  for  the  stress  assisted  diffusion  falls  as  KT  falls.  Both 
effects  decrease  the  amount  of  hydrogen  reaching  the  plastic  zone  and 
hence  decrease  the  rate  of  hydride  formation  resulting  in  a  decrease  of 
crack  velocity  as  observed  in  Region  C  (Fig.  60c).  The  crack 
propagates  at  a  decreasing  velocity  until  the  stress  intensity  falls 
below  K^q>  the  stress  intensity  to  crack  the  hydride,  or  until  Kj 
falls  below  the  stress  required  to  form  new  hydrides  at  the  temperature 
of  the  test.  The  latter  factor  would  provide  the  threshold  condi¬ 
tion  at  temperatures  above  the  solvus  while  the  former  would  apply 
below  the  solvus  since  even  relatively  small  stress  intensities  would 
bias  the  free  energies  in  favor  of  hydride  formation  at  the  crack 
tip.  This  expectation  is  in  accord  with  the  data  presented  in  Figs. 
31-34.  The  data  for  d.t.  specimens  23  and  25  (Fig.  33)  correspond  to 

the  tests  95  K  below  the  solvus  and  the  in  Region  C  is  about 
1/2 

3  MPa  •  m  .In  contrast,  the  tests  for  d.t.  specimens  53  (Fig.  31), 

41  and  54  (Fig.  34),  which  correspond  to  test  temperatures  13,  18  and  16 

K  above  their  respective  solvus  temperatures,  have  threshold  stress 

1/2 

intensities  of  about  7,  8  and  11  MPa  •  ra  ,  respectively.  These 
would  correspond  to  the  required  to  shift  the  local  solvus  tempera¬ 
tures  to  the  test  temperatures  and  are  larger  than  as  expected. 


Further  discussion  of  K_  for  various  tests  will  be  found  in  section 

Io 

6.4 

The  data  exhibiting  three  region  behavior  discussed  above 
was  obtained  in  decreasing  tests  at  intermediate  temperatures.  A 
similar  behavior  was  noted  in  the  decreasing  K^  test  at  112  K  (Fig.  25) 
although  the  data  were  not  as  complete.  This  test  was  performed  96  K 
below  the  solvus.  While  the  general  behavior  was  similar  to  chat 
observed  for  specimens  tested  in  the  decreasing  test  at  148  K,  95  K 
below  the  solvus  (Fig.  33),  the  three  regions  were  less  well 
defined.  Double  torsion  specimens  tested  at  148  K  had  almost  20  times 
more  hydrogen  in  solution  than  in  the  specimen  tested  at  112  K  (see 
Table  6).  It  appears,  however,  that  the  small  amount  of  hydrogen 
remaining  in  solid  solution  for  the  specimen  tested  at  112  K  is 
sufficient  to  cause  the  three  region  (V-K)  curve  though  to  a  much  less 
well  defined  extent. 

The  distinguishing  feature  of  the  intermediate  temperature 
range  fractography  is  the  completely  brittle  fracture  in  Regions  A 
and  C  and  the  large  amount  of  ductility  associated  with  Region  3  o' 
the  (V-K)  curves  (Figs.  39-43).  This  is  consistent  with  the  suggestion 
that  the  increase  in  crack  velocity  for  decreasing  in  Region  3 
is  due  to  the  hydride  formation  in  front  of  the  crack  and  linking  ut> 


of  fracturing  hydrides.  A  large  amount  of  ductility  between  these 
hydrides  would  be  expected  if  the  velocity  was  too  large  for  the 
hydrogen  flux  which  results  and  therefore  a  complete  hydride  path  is 
not  formed.  The  more  brittle  f ractographic  features  for  Regions  A 
and  C  as  shown  previously  (Figs.  36-38  and  44)  are  consistent  with 


their  corresponding  (V-K.)  behavior,  in  Region  A,  the  precipitation  of 
the  hydrides  during  loading  results  in  large  hydrides  which  form  an 
almost  continuous  path  for  the  fracture.  In  Region  C,  the  rate  of 
fracture  is  low  enough  and  is  controlled  by  hydride  formation  at  the 
crack  tip  at  a  rate  sufficient  to  form  a  complete  hydride  path  for 
the  fracture. 

6.2.2  Constant  Data 


The  occurrence  of  the  three  stage  behavior  for  the  constant 
tests  as  shown  in  Fig.  45  can  be  readily  explained  by  a  stress 
induced  hydride  embrittlement  mechanism.  In  the  low  Stage  I  an 
extremely  high  dependence  of  crack  velocity  on  K^.  was  observed. 

The  nature  of  this  dependence  is  such  that  Stage  I  may  be  a  transi¬ 
tion  between  the  crack  arrest  KT  (i.e.  K_  ,  the  stress  intensity  for 

l  Io 

crack  propagation  in  the  brittle  hydride)  and  the  constant  ,  Stage 

II  behavior.  This  is  in  agreement  with  the  work  done  on  hydrogen 

(27) 

embrittled  Zirconium  alloys.  Discussion  of  K_  for  constant  KT 

Io  I 

will  be  found  in  section  6.4 


Stage  III  was  characterized  by  the  onset  of  rapid  crack  growth 

and  was  observed  to  be  independent  of  temperature  and  hydrogen 

concentration.  It  corresponds  to  the  effective  critical  stress 
0  f  f 

intensity,  ,  for  rapid  fracture  by  overload.  The  effective  critical 

stress  intensity  will  be  less  than  the  true  since  the  presence  of 

birttle  hydrides  in  a  ductile  matrix  will  reduce  the  fracture 


roughness.  Fractography  was  consistent  with  this  behavior  since  the 
amount  of  ductility  increased  as  Stage  III  was  approached  (Fig.  50). 
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Stage  II  for  the  constant  data  can  be  explained  by  a  rate 

limiting  model  having  the  crack  velocity  controlled  by  the  rate  of 

diffusion  of  hydrogen  to  the  crack  tip.  In  this  region  a  continuous 

hydride  fracture  path  forms  at  a  rate  which  is  controlled  by  hydrogen 

flux  to  the  crack  tip.  A  modified  theoretical  model  by  Simpson 
(27) 

and  Puls  based  on  diffusion  limited  crack  propagation  is  presented 

in  section  6.3.  This  model  allows  the  expected  activation  energy  for  the 
rate  limiting  process  of  Stage  II  to  be  calculated.  The  theoretical 
value  will  be  shown  to  be  in  good  agreement  with  the  experimentally 
observed  activation  enthalpy  for  crack  propagation  in  Stage  II  of  the 
constant  data. 


6 . 3  Theoretical  Stage  II  Activation  Energy 


(27) 

A  theoretical  model  was  presented  by  Simpson  and  Puls  to 

describe  the  crack  propagation  behavior  in  hvdrogen  embrittled  Zr-2.5 

wt.%  Nb  alloys.  A  cylindrical  geometry  shown  in  Fig.  61  is  assumed 

( 26 ) 

surrounding  the  crack  tip.  Their  expression  for  the  crack  velocity 

when  applied  to  Nb-H  alloys  is  given  by 


V 
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where 


W 


a 

t 


^hyd’^B 


No 

H 

o 

H 


L 


l 


V 

*(i,l) 


Total  molal  elastic  strain  energy  of  matrix  and 
hydride. 

Molal  interaction  energy  of  an  inclusion  with  an 
applied  elastic  stress  field. 

Atomic  volumes  of  hydride  and  alpha  phase  atoms  in 
Nb ,  respectively. 

Diffusion  coefficient  of  hydrogen  in  alpha  Mb. 

Terminal  solid  solubility  of  hydrogen  at  zero 
applied  stress  when  precipitated  hydrides  are  present 
(or  the  average  hydrogen  concentration  in  solid 
solution  at  temperatures  above  the  solvus). 

Distance  from  the  crack  tip  to  source  hydrides  in 
Zr-H  alloys  as  shown  in  Fig.  61.  This  parameter  was 
redefined  below  because  source  hydrides  do  not 
exist  in  Nb-H  alloys. 

Distance  from  crack  tip  to  center  of  plastic  zone  given 
7  2 

by  K“/l2TTOy  .  The  factor,  a  ,  is  the  yield  scress 
and  at  low  temperatures  is  a  strong  function  of 
temperature . 

Crack  tip  radius  and  is  given  by  K^/2a.E,  where,  E, 
is  Young's  Modulus. 

Function  of  l  and  L  approximated  by  ln(L/2). 


Equation  (34)  was  derived  assuming  that  the  rate  of  hydride  growth  in 
the  vicinity  of  the  crack  tip  is  equal  to  the  crack  velocity.  This 
ignores  the  observation  chat  the  hydrogen  embrittled  crack  growth  is 


61.  Geometry  of  crack  and  hydrides  used  in  deriving  Eqn.  (34). 
(After  R.  Dutton  et  al,).(26) 


119 


generally  discontinuous  in  nature.  Since  resolution  of  the 

hydrides  in  niobium  are  not  observed  as  the  crack  propagates,  the 
concept  of  hydrides  acting  as  sources  of  hydrogen  does  not  apply 
and  the  parameter  L  must  be  redefined.  In  the  Nb-H  system  the  solid 
solution  acts  as  a  source  of  hydrogen  for  the  precipitation  of  the 
hydrides  at  the  craufc.  tip.  In  this  case  the  parameter  L  is  the 
distance  at  which  the  hydrogen  concentration  equals  the  solid 
solution  concentration  in  equilibrium  with  the  hydrides  in  an  unstressed 
system  (or  the  solid  solution  concentration  in  the  absence  of  the 

_3 

hydrides).  For  computational  purposes,  L  was  taken  to  be  5  x  10  m 

which  is  approximately  ten  times  the  plastic  zone  radius.  The  crack 

velocity  depends  only  on  [ln(L)]  1  since  W^nc  does  not  depend 

sensitively  on  the  parameter  L  (Eqn.  (34)). 

The  major  temperature  dependence  in  Eqn.  (34),  neglecting 

W*nc  and  Wa  is  found  in  the  product  D^3  •  C°  for  specimens  tested 
t  t  Hr! 

below  the  solvus  temperature  in  systems  where  the  hydrides  can  redissolve 

and  act  as  a  source  of  hydrogen.  For  tests  conducted  above  the 

solvus  temperature  or  where  resolution  of  the  hydrides  does  not  occur, 

Mb 

temperature  dependence  resides  only  in  D  .  The  diffusion  coefficient 

Cl 

can  be  expressed  as 

Mb 

Dh  =  D^exp (-Q/RT)  ,  05) 

where  «  the  pre-exponential  multiplier, 


Q  =  the  activation  energy. 
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Neglecting  W^nc  and  Eqn.  (34)  can  be  rewritten  as, 

V  =  aC°exp(-Q/RT)  (36) 

where 


»  *  (451hyd^Do)/3(J,')aNb$(2,,L)  .  (37) 


The  crack  velocity  given  in  Eqn.  (36)  should  follow  an  Arrhenius 
behavior.  The  diffusion  coefficient  for  hydrogen  in  niobium,  measured 

(45) 

by  Schaumann  et  al.  has  been  found  to  deviate  from  an  exponential 

behavior  below  room  temperature.  On  the  other  hand,  Matusiewicz  et 
(41) 

al.  did  not  find  this  deviation  and  for  temperatures  greater  than 


165  K,  the  diffusion  coefficient  had  an  activation  energy  of  8.20  kJ/mol. 

The  factor  W*nc  depends  only  on  the  degree  of  elastic  misfit 
between  the  matrix  and  hydride  but  not  on  the  applied  stress.  This 


implies 


(27) 


,.inc  / « \  -.itiCxx  .  _xr\c 
W  (Z)  =  Wt  (L)  * 


which  allows  its  contribution  to  the  activation  energy  to  be 

approximated.  Puls^  calculated  W^nc  for  delta  and  gamma  Zr  hydrides 

having  various  orientations.  Values  he  found  for  hydrides  having 

shapes  of  plates  and  needles  with  various  orientations  ranged  f rom  7.74 

kJ/mol  to  2.66  kJ/tnol.  Morphologies  observed  for  niobium  hydrides  were 

(2) 

plate,  blocky,  dendritic  and  starlike.  Stress  induced  hydrides 

appearing  on  specimen  surfaces  resembling  ellipsoids  were  observed  in 
this  investigation.  Knowing  the  equilibrium  shape  of  niobium 
hydride  i.e.  a  plate,  lying  parallel  to  its  { 100 }  fracture  plane, 
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for  niobium-hydrogen  alloys  was  calculated  by  the  following 
expression  used  by  Puls. ^ 

Winc  -  Vq[E/(1-v)]{A2  +  (fTc/2a)  [A*  £  +  £2/4  (1+v)  ] }  (39) 

t  p 

where  c/a  =  0.05  (ratio  between  minor  and  major  axes  for  plate). 

E  *  1.05  x  10  ^  Pa  (Young's  modulus  for  niobium)  . 

\i  *  0.39  (Poisson's  ratio  for  niobium)/4^ 

—6  3 

*  12.3  x  10”  mJ/mol  (partial  molal  volume  of  beta 
p 

hydride  calculated  from  lattice  constants  for  niobium 
hydride) . 

A  and  £  *  Values  defined  by  the  transformation  (stress  free)  strains 

of  beta  hydride  which  were  calculated  from  lattice 
(47) 

constants  and  are  equal  to  0.039  and  0.003, 
respectively. 

Substituting  these  values  into  Eqn.  (39)  yields 

Winc  *  3.24  kJ/mol  (40) 

t 

i.nc 

Since  W  is  positive,  increasing  values  will  exponentially  increase 
the  crack  velocity. 

The  factor  W*  in  Eqn.  (34)  can  be  written,  assuming  elastic 

(5) 

constants  of  hydride  and  matrix  are  equal,  as 


W 


a 

t 


'Vijeij 


(41) 


where 


applied  stresses  (positive  when  tensile), 

stress  free  transformation  strains  for  the  alpha  to 


beta  transition. 


It  can  be  shown  that 


W®(1)  »  U*(L)  (42) 

since  the  hydrostatic  stresses  at  the  crack  tip  are  much  greater  than 
the  stresses  outside  the  plastic  zone.  Typical  stresses  occurring  at 
l  (plastic  zone)  under  plane  strain  conditions^  are 


a 


11 


3oy(T)  , 


°22  “  1-8°y(T)  * 

°33  *  2.4ay(T)  . 


(43) 


Using  the  data  from  Chen  and  Arsenault^*^  and  Pereturina  and  Pavlov,^ 
the  temperature  dependent  yield  stress  of  hydrogen  charged  polv- 
crystalline  niobium  was  approximated  by 


a  (T)  -  (1.09  x  105  MPa  *  K_1)/T  -  (T/4)MPa  •  K"1  ,  (44^ 

for  temperatures  between  75  and  300  K. 

Figure  62  is  W*(2.)  versus  temperature  for  a  hydride  plate 
parallel  to  {100}  and  the  stresses  and  applied  along 

[001],  [100]  and  [010]  directions,  respectively.  For  temperatures 
between  180  and  303  K,  which  correspond  co  the  constant  data, 

W^( 1)  ranges  between  -2.0  and  -1.0  kJ/raol. 

Including  W*(£)  and  W^nc ,  Eqn.  (36)  becomes 


(kJ/mol) 
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Temperature  (K) 

Fig.  62.  Graph  of  -W®  versus  temperature  for  niobium  hydride  plate 
parallel  to  the  (100)  planes  and  the  stresses  given  by 
Eqn.  (43)  applied  along  the  [001],  [100]  and  [010]  directions, 
respectively. 
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Substituting  the  values  for  Q,  4Tnc  and  W*(Jl)  into  Eqn.  (45),  the 
total  activation  energy,  Q",  is  found  to  range  between 

5.96  to  6.96  kJ/tnol  (46) 

for  the  constant  data.  This  range  of  activation  energy  is  within 
the  experimental  error  of  the  measured  temperature  dependence  as  will 
be  shown  below. 

The  Stage  II  activation  enthalpy  of  the  experimental  constant 
iCj  data  can  be  compared  to  the  theoretical  values  of  Eqn.  (45).  Before 
this  can  be  done,  the  effect  of  the  hydrogen  concentration  in  solution 
on  the  crack  velocity  must  be  taken  into  account.  For  the  Nb-H  system 
the  appropriate  value  of  C°  is  the  concentration  in  solid  solution. 

rl 

This  is  achieved  by  rewriting  Eqn.  (45)  so  chat 

V/C°  -  aexp(-AH/RT)  (47) 

where  AH  is  the  activation  enthalpy  of  the  Stage  II  constant  data. 
Taking  the  natural  logarithm  of  Eqn.  (47)  yields 

ln(V/C°)  -  lnct  -  AH/RT  .  (48) 

n 

The  BLj.  range  common  to  Stage  II  behavior  for  all  temperatures 

1/2 

tested  lies  between  27.5  and  32.5  MPa  •  m  (Fig.  45).  Figure  63 

1  > 2 

shows  the  Arrhenius  plot  of  Eqn.  (48)  for  *  30.0  MPa  •  m  . 

Similar  behavior  was  observed  for  ■  27.5,  32.5  and  35  MPa  •  “. 

Using  Che  following  relation, 

AH  -  R[Aln(V/C°)/A(l/T)]  (49) 


-ln(V/C 


11 


T 


12 


Kx  =  30  MPa/rrT 


1g  I _ i _ I _ i _ I _ i _ L 

3  4  5  6 

1000/T  (K'1) 


Fig.  63.  °lot  of  -In(V/C^)  versus  l'T  for  -  30  MPa  •  rr.‘  . 

activation  enthalpy  from  best  fit  to  data  (solid  ;:nv 
Eqn.  ( 49 )  is  6.65  k.J  'mol. 


where  R  =  the  gas  constant, 


T  »  temperature  in  K, 

the  activation  enthalpies  were  calculated  and  then  averaged  resulting 
in 


2.H  *  6 . 6  r  1 .  5  kJ/mol  (50) 

which  agrees  very  well  with  the  measured  values  given  in  Eqn.  <  36)  . 

Figure  b~>  shows  Eqn.  (3d)  plotted  as  crack  velocitv  versus  R. 
for  temperatures  of  200,  230  and  J00  K.  Also  plotted  for  purposes  it 
comparison,  are  the  constant  Kr  curves  tor  temperatures  of  2  32  ano  30 
The  factor  1‘  appearing  in  the  denominator  of  F.qn.  .  3d  i  was  defined 

~i 

earlier  as  being  proportional  to  Rjl  •  As  seen  in  Fig.  h-  ,  the  1 V — K  ' 

curves  of  Eqn.  (  3*  i  reflect  the  crack  velocities  major  dependence 

(l'K^i.  The  functional  dependence  if  Eqn.  (34)  does  not  resemble  tve 

i  1 7  ' 

experimental  lv  observed  curves  for  either  T.r-H  "  or  Xh-H.  illoo- 
>Fig.  n- '  for  the  -nst^nt  Kr  Stige  11  oehuvior  and  require,  •  .r'n-.  r 
refinement.  Numer  »  ,;a  1 1  v  ,  nowever.tne  cal  .1  it  eg  ve..  ••  :  t  t«-s  ir  w  : 
i  fact  >r  of  three  ■  :  toe  St  ige  11  .oast ant  ■.  -racs  '.el.  it  .es  :  : 

equal  R  values. 

n. «  Th res no  Id  Stress  Intensttv  Fa-  tors 
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stresses  on  the  solvus  as  was  briefly  reviewed  in  the  Previous 
Research  section.  For  specimens  tested  below  its  solvus,  hydrides 
were  present  before  loading  occurred  and  the  >  Tg)  is  related 

to  the  minimum  load  required  to  fracture  the  existing  hydrides  and 
form  stress  concentrations  at  which  additional  hydrides  are  stress 
induced.  Repetition  of  this  process  results  in  slow  crack  growth. 
Specimens  tested  above  their  solvus  temperatures  did  not  have  hydrides 
present  prior  to  loading  and  as  discussed  in  section  6.2.1  for  the 
decreasing  K7  tests,  crack  propagation  will  not  occur  until  the 
stress  intensity  and  the  crack  tip  plasticity  raises  the  local  solvus 
temperature  to  the  test  temperature  and  allows  hydride  precipitation. 
The  difference  between  the  stress  intensitv  thresholds  above  and  below 
the  solvus  temperatures  represents  the  additional  work  required  to 
stress  induce  hvdrides.  Similar  behavior  for  the  threshold  stress 
intensities  was  observed  for  the  decreasing  kT  tests  as  shown  in 

L 

Figs.  31-35  which  show  K„  1 T  T  )  greater  than  K,  (T  <  T  ) . 

io  s  *  Io  s 

The  constant  KT  specimen  tested  at  the  highest  temperature. 


i.e.,  303  K,  was  1-*  K  above  its  solvus  and  had  a  slightly  larger  stress 

intensitv  threshold  than  trie  other  two  specimens  tested  jC  10  k  and  13 

above  their  solvus  temperatures.  This  is  consistent  with  the  above 

model  since  tie  h liner  i  specimen  is  tested  above  Its  solvus  t-mper iture 

tie  greater  the  threshold  stress  intensitv  should  be. 

Fi  g  :re  ;  s  i  gran'.  :  K.  '  T  ■  -ersis  temper  at  are  for 

i  ■>  s 

.n  r»-  ising,  ie  Teasing  and  ns:  ml  K,  tests.  The  increasing  me 
•••  -  •  i  -  ’  *.  "-res-:  in  - «  ep.  t  •  •  •  1 1  ow  the  same  trend,  i.e.. 

•  .  ’  .  n  re  is  es  w.  •  •  r-'ismg  temper  it  are.  The  single  point 


Threshold  Stress  Intensity  Factor  (MPa 
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Temperature  (K) 


65.  Threshold  stress  intensity  factor  below  the  solvus  versos 
temperature  for  increasing  (circles),  constant 
(square)  and  decreasing  K^.  (triangle)  tests.  Solid  line  is 
least  squares  fit  to  increasing  K  and  constant  K  data. 
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plotted  for  the  decreasing  tests  below  its  solvus  indicates  a 

1/2 

K  (T  <  Tg)  of  only  3  MPa  •  m  “at  148  K  and  differs  significantly 

from  the  increasing  and  constant  K  tests. 

Gahr  et  al. calculated  the  surface  energies  for  cleavage 

of  niobium  hydride  single  crystals  ac  77  K  and  found  the  critical  stress 

1/2 

intensities,  K^,  Co  range  rrom  0.9  to  1.1  MPa  •  m  .  This  indicates 

that  additional  factors  are  responsible  for  the  larger  K.T  (T  <  T  ) 

lo  s 

observed  for  the  various  tests  in  the  present  studies.  The  hydrides 
precipitated  on  cooling  below  the  solvus  in  the  fractured  specimens 
were  surrounded  by  a  ductile  solid  solution  matrix.  Due  to  the 
volume  increase  on  forming  the  hydride,  the  precipitates  are  in 
compression  while  the  surrounding  solid  solution  is  highly  deformed . ^ ^ 
Under  these  circumstances  the  K^fT  <  Tg)  may  reflect  the  stress 
intensity  to  form  dislocation  pile-ups  in  the  solid  solution  which 
initiate  the  cracking  in  the  hydrides.  Thus  in  constant  KT  and 
increasing  KT  tests  the  K_  (T  <  T  )  reflects  the  plastic  oroDerties  of 
the  solid  solution  surrounding  the  precipitated  hydrides.  As  a 
result  Kjq(T  <  Tg)  is  larger  than  the  K^o  measured  on  single  crystals 
of  hydrides  and  exhibits  a  strong  temperature  dependence  (Fig.  65). 

In  contrast  to  this  behavior  the  low  value  observed  (Fig.  65)  in 
decreasing  tests  corresponds  to  the  hydride  since  in  this  test 

the  K^o  measures  the  stopping  of  a  crack  which  is  already  propacating  in 
a  hydride  formed  at  higher  Kj.  values.  Since  the  crack  aireadv  exists 
in  the  hvdride  the  stress  required  to  stop  its  propagation  in  the 
decreasing  KT  test  should  correspond  to  the  value  measured  bv  ,'uh r 


eC  al.  The  value  of  about  3  MPa  •  ra^“  is  in  relatively  good 

1/2 

agreement  with  the  value  of  1.0  MPa  •  m  measured  for  the  single 


crystal  niobium  hydride. 
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7.  SUMMARY  AND  CONCLUSIONS 

The  underlying  mechanism  for  the  observed  fracture  behavior 
of  hydrogen  embrittled  niobium  was  shown  to  be  consistent  with  the 
fracturing  of  preexisting  and  stress  induced  hydrides.  The  specific 
behavior  of  the  various  (V-K)  curves  was  shown  to  depend  on  the  type 
of  loading  a  specimen  underwent.  Similarities  and  differences  in 
the  (V-K)  behaviorwere  observed  for  increasing,  decreasing  and  constant 
K-j.  tests. 

Results  of  the  decreasing  K^  tests  were  able  to  be  grouped  in 
low,  high  and  intermediate  temperature  ranges.  The  relationship 
between  the  specimen's  test  and  solvus  temperature,  proved  to  be 
the  determining  factor  in  the  type  of  (V-K)  behavior  observed.  A 
unique  three  region  (V-K)  behavior  occurred  in  the  intermediate 
temperature  tests.  These  regions  were  denoted  A,  3  and  C  and  are 
distinctly  different  from  the  three  stage  behavior  observed  in  tne 
constant  Kj  tests.  The  unusual  observation  of  the  crack  velocity 
increasing  with  decreasing  K^  (Region  B)  was  shown  to  be  due  to  the 
linking  up  of  fractured  hydrides  to  form  one  large  crack. 

Differences  in  the  stress  intensity  threshold  for  the  constant 
Kj  specimens  tested  above  and  below  their  solvus  temperatures  were 
successfully  explained  in  terms  of  the  additional  stress  required  to 
raise  the  solvus  tc  allow  hydride  precipitation. 

Comparison  between  the  experimental  and  theoretical  activation 
energies  for  Stage  II  of  the  constant  K^  tests  showed  excellent 
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agreement.  This  gives  further  support  to  the  proposed  model  of  hydrogen 
diffusion  to  the  crack  tip  being  the  rate  limiting  mechanism  for  crack 
growth  during  Stage  II.  Comparison  between  crack  velocities  measured 
in  hydrogen  and  deuterium  charged  specimens  used  in  decreasing  K^  tests 
also  showed  excellent  agreement  with  the  expected  differences  due  to 
the  isotope  effect.  This  again  indicated  that  crack  velocities  are 
related  to  the  dif f usivities  of  hydrogen  or  its  isotope. 

Fractographv  was  able  to  be  correlated  with  (V-K)  behavior 

of  the  various  tests.  Fracture  surfaces  of  the  s.e.n.  specimens  showed 

(2) 

identical  features  of  previous  work  where  transgranular  cleavage 
was  the  only  brittle  failure  mode  observed.  Fracture  surfaces  of 
intermediate  temperature,  decreasing  Kj  tests  reflected  its  corresponding 
three  region  (V-K)  curves  particularly  well.  Brittle  behavior 
occurred  at  low  crack  velocities  and  more  ductile  features  were  present 
at  higher  velocities  reflecting  the  linking  up  of  fracturing  hydrides. 


APPENDIX  A 


PHOTOELASTIC  INVESTIGATION  OF  DOUBLE  TORSION  GEOMETRY 

A. 1  Introduction 

Photoelastic  investigation  of  the  d.t.  specimen  was  undertaken 

in  order  to  document  the  stress  distribution  in  the  plane  of  the 

specimen.  No  analysis  of  this  stress  distribution  appears  in  the 

(33) 

literature  although  Trantina  used  a  three  dimensional  elastic 
finite  element  stress  analysis  to  measure  in  the  plane  of  the 
crack.  For  models  which  are  loaded  to  produce  generalized  plane 
stress  conditions,  the  two  principal  stress  directions  and  magnitudes 
do  not  change  appreciably  with  model  thickness.  As  a  result,  classical 
photoelastic  techniques  could  be  used  for  their  stress  analysis. 

Though  the  d.t.  specimen  is  in  the  shape  of  a  rectangular 
plate,  a  three  dimensional  photoelastic  technique  v.i.z.  the  "frozen 
stress"  method  was  used  due  to  its  particular  loading  arrangement. 

This  method  is  normally  reserved  for  more  complicated  configurations. 
Stressed  models  whose  three  principal  stress  axes  change  in  magnitude 
and/or  inclination  in  such  a  way  that  they  are  parallel  to  the  incident 
ray  of  polarized  light,  require  the  "frozen  stress"  technique.  Bv 
using  the  proper  polymer  and  load-thermal  cycle,  the  fringe  pattern 
from  an  applied  stress  can  be  "frozen"  into  tne  model.  The  model 
can  subsequently  be  cut  into  thin  slices,  without  altering  the  fringe 
pattern,  and  examined  as  a  generalized  plane  stress  model. 
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As  discussed  by  Hendry^^  the  phenomenon  of  stress  freezing 
is  due  to  the  existence  of  two  phases  in  the  plastic.  One  is  very 
temperature  dependent  being  rigid  at  room  temperature  and  fluid  at 
high  temperature.  The  other  phase  has  a  small  temperature  dependency  and 
is  rigid  at  all  temperatures  but  less  so  than  the  low  temperature 
rigidity  of  the  first  phase.  A  mechanical  model  consisting  of  a  spring 
in  a  fluid  which  is  rigid  at  room  temperature  but  fluid  at  the  softening 
temperature  represents  the  stress  freezing  behavior  quite  well. 


A. 2  Experimental  Procedure. 

The  plastic  d.t.  specimen  was  machined  out  of  a  1.02  cm  thick 
plate  of  precast  epoxy  resin  called  PSM-9  obtained  from  Photoelastic 
Division,  Measurements  Group  in  Raleigh,  North  Carolina.  The  specimen 
was  5.6  cm  wide  and  12  cm  long.  A  notch  4  cm  long  on  top  and  3.2  cm 
on  bottom  of  the  specimen  was  cut  with  a  jeweler's  saw.  After  machining, 
the  specimen  was  placed  in  the  Instron  loading  stage  and  immersed  in 
an  inert  silicone  oil  bath.  Heating  and  cooling  rates  were  controlled  by 
passing  the  current  from  a  clock  driven  Variac  through  a  400  W  heater. 
Temperature  was  monitored  by  thermocouples  and  kept  uniform  by 
circulating  the  oil.  The  heating  rate  was  150  K/hr.  At  a  final  temper¬ 
ature  of  130  K,  a  load  of  10  N  was  applied  to  the  specimen.  A  crack 
initiated  and  advanced  approximately  3  cm  before  arresting.  While 
maintaining  this  load,  the  specimen  was  annealed  for  2  hours.  Tempera¬ 
ture  was  then  lowered  to  room  temperature  at  a  rate  of  2  K/hr.  After 
the  stress  freezing  cycle,  the  specimen  was  sliced  in  the  manner  shown 
in  Fig.  66.  The  slices  consisted  of  six  squares  3.3  cm  on  a  side 


corresponding  to  six  tenths  of  the  specimen  width.  Thicknesses  are 
listed  in  Table  9.  Slices  two  and  three  each  consist  of  two  halves  since 


the  propagating  crack  bisected  them.  Slices  one,  five  and  six  were 
partially  bisected.  Slices  one  and  four  represent  the  top  surface 
of  the  specimen,  two  and  five,  the  middle  and  three  and  six,  the  bottom. 
These  slices  were  cut  using  a  water  cooled  jewelers  saw  to  minimize 
generation  of  heat  which  would  tend  to  distort  the  fringe  pattern  at  the 
cut  edges.  Cut  surfaces  were  wet  sanded  to  eliminate  surface 
scratches.  A  circular  polariscope,  model  series  241  manufactured  by 
the  supplier  of  PSM-9,  was  used  to  observe  the  isochromatics  and 
isoclinics.  A  fringe  multiplier  was  used  to  obtain  reasonable  numbers 
of  observable  isochromatics .  Only  a  multiplication  of  five,  however, 
was  attainable  before  off  axis  distortion  prevented  the  resolution 
of  higher  order  fringes.  All  isochromatics  were  photographed  in 
bright  field  on  Kodak's  Contrast  Process  Ortho  4154  film  using 
monochromatic  light  from  a  green  filtered  mercury  lamp.  Isoclinics 
were  photographed  on  PN  55  Polaroid  film  using  a  diffuse  fluorescent 
light  source.  Isoclinics  for  the  middle  section  were  difficult  to 
locate  and  their  presence  for  angles  between  40°  and  90°  had  to  be 
estimated.  The  middle  section  only  exhibited  two  orders  of  isochro¬ 
matics  even  under  a  fringe  multiplication  of  five. 

A . 3 .  Model  Fringe  Value  Calculation 

A  measurement  of  the  material  stress  fringe,  f,  for  PSM-9  was 
necessary  before  the  principal  stress  differences  could  be  calculated 
from  the  isochromatics.  The  fringe  value  is  defined  as  the  stress 


U 


required  Co  produce  n  number  of  fringes  in  a  place  of  thickness  t  jnj  i, 

,  (52) 
given  by 


f  *  ( j  -  ^2) c/n  ( 31  1 

A  calibration  disk,  loaded  in  compression  along  a  diameter  was 

used  co  find  the  f  value.  The  disk  was  cut  to  a  diameter  of  3.39  -m 

on  a  precision  lathe  out  of  the  same  sheet  used  in  making  the  oiasti.. 

d.t.  model.  The  thermal  cycle  of  the  disk  was  controlled  tc  closelv 

duplicate  the  thermal  cycle  of  the  d.t.  model.  Isoehromacics  of  the  disk 

were  photographed,  without  a  fringe  multiplier,  in  green  monochromatic 

light.  Both  bright  and  dark  field  were  used  in  order  to  obtain  one 

half  fringe  resolution.  The  analytic  expression  for  the  principal 

(53) 

stress  differences  at  the  center  of  the  disk  is  given  bv' 

~L  -  C,  =  8P/Tdt  (32) 

where  P  *  applied  compressive  load 
d  =  diameter  of  disk 

At  the  center  of  the  disk,  the  fringe  order  was  found  to  be  -*.8 
for  an  applied  load  of  32.5  N.  Substituting  Eqn.  (52)  into  Eqn.  (51) 
and  solving  for  f  yields  f  ■  3.20  x  10  “*  Pa'fr/cm.  This  value  must  then 
be  divided  by  five  to  obtain  the  proper  material  fringe  value  for  the 
slices  since  their  isochromatics  were  obtained  using  a  fifth  order  ‘‘ringe 
multiplication,  so 

f  ..  ■  6. -*0x10^  Pa/fr/cm  . 

slice 


(5  3) 


The  a*  ..lei  :  r :  rjge  i :  ue  i.  r  -r  >  . .  i  «  ir«-  .  utfo  .• 

Tinie  ■/ . 

A.-*  Principal  Stress  Ji  i  ferenccs  u-.g  ':rc  t 

In  this  section,  tne  i  sncn  r  mut  i  :s  me  iso..  ;  m  i  s  :  r  t  .•  • 

miudie  and  bottom  slices  ire  given.  The  principal  ->:res-  :  .  :  •  ;  ns 

ire  drawn  f  :  r  the  isoelinios  'using  trie  me  mods  :  .'ess.  0  no  it.:  r  r  1  -  , 

mil  the  principal  stress  differences  for  ir'.vri.  sri,-  *eo  dir.  ns 

through  the  slices  ire  compared.  It  is  important  to  :n  te  mat  : 
calculated  principal  stresses  ire  not  the  nr  in-  ipa.  stresses  hut 
instead  projections  of  the  principal  stress  .1  i :  :  erences  in  re 

the  slices.  With  the  information  given  below,  it  is  not  rioss;:  .e  t 
resolve  the  stress  differences  into  their  principal  componer.  ts .  A 
comprehensive  solution  to  the  stress  distribution  in  the  d.t.  iniv  iten 
would  require  1  second  identically  prepared  plastic  d.t.  specimen  mode 
wnose  slices  were  perpendicular  to  the  original  six.  Then,  in 
theory,  the  principal  stresses  could  be  resolved  into  their  in.i .  v  1  iu.i  1 
components.  In  practice,  this  is  j  very  tedious  procedure.  The 
following  analysis  is  not  intended  for  exact  work  but  is  presented  t 
an  approximate  indication  of  the  relative  magnitude  of  the  projected 
principal  stress  differences  in  the  top,  middle  and  bottom  slices. 

Figures  67-69  show  the  isoelinios  and  stress  tr  1  i ee t >r  ies .  ii 
the  slices  as  well  as  the  entire  specimen  are  symmetrical  about  the 
plane  of  the  crack  surface,  the  isoclinics  are  drawn  on  the  left  half 
and  the  principal  stresses  are  drawn  on  the  right  half  for  claritv  of 
presentation.  Figures  70-72  show  the  isochromatics  tor  the  ton,  midi 


143 


Fig.  68.  Isoclinics  (dashed  curves)  and  principal  stresses  (solid 
curves)  for  middle  slices  (two  and  five). 


Fig.  69.  Isoclinics  (dashed  curves)  and  principal  stresses  (solid 
curves)  for  bottom  slices  (three  and  six) . 
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Fig.  70.  Isochromatic  fringes  for  top  slices  (one  and  four). 

Lines  0A,  OB  and  0C  are  graduated  in  tenths  of  (W/2) 
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5/2 
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and  bottom  slices,  respectively.  All  curves  represent  odd  half  numbers 
of  relative  wavelength  retardation  and  are  labled  for  fringe  order 
where  space  allowed.  Also  drawn  on  the  right  half  of  these  isochromatic 
figures  are  three  parallel  lines  labeled  OA,  OB  and  OC  which  lie 
perpendicular  to  the  crack  direction.  Each  line  is  graduated  in 
fractional  units  of  one  half  the  width,  W  of  the  specimen.  OA  inter¬ 
sects  the  crack  tip  on  the  bottom  slice,  OC  intersects  the  crack  origin 
on  the  top  slice  and  OB  intersects  the  fracture  surface  on  Che  top, 
bottom  and  middle.  It  is  of  interest  to  examine  how  the  magnitude  of 
the  principal  stress  difference  varies  along  lines  OA,  OB  and  OC 
for  the  three  different  levels.  These  results  are  shown  in  Figs. 

73-75  and  are  discussed  below. 

The  approximate  state  of  stress  near  the  crack  front  can  be 
recognized  simply  by  considering  the  loading  arrangement  of  the  d.t. 
specimen,  thus,  the  crack  front  appearing  on  the  top  surface  should  be 
in  a  state  of  compression.  The  crack  surface  in  the  middle  section 
should  be  in  a  type  of  neutral  stress  zone,  similar  to  the  neutral 
plane  which  occurs  in  a  beam  loaded  in  pure  bending.  On  the  bottom 
of  the  specimen,  the  crack  surface  should  be  in  a  state  of  tension. 

Referring  to  Fig.  73,  the  principal  stress  difference  along 
OA  intersects  the  crack  tip  in  the  bottom  surface  and  therefore 
exhibits  the  greatest  stress  concentration.  The  fractional  orders 
of  the  isochromatics  intersecting  OA  for  the  middle  section  were  too 
ill  defined  to  locate,  however,  they  were  substantially  less  chan  the 
stresses  for  the  top  or  bottom.  At  the  leading  edge  of  the  crack 
front,  which  corresponds  to  the  origin  of  line  OA  on  the  bottom 


16 


Fig.  73.  Principal  stress  difference  versus  fractional  distance  along 
line  OA  for  top  and  bottom  slices. 


Stress  (10'^  Pa) 


0  0.1  0.2  0.3  0.4  0.5  6.0 

Oistance  X-W/2  Along  Line  OB 


Fig.  74.  Principal  stress  difference  versus  fractional  distance  along 
line  OB  for  top,  middle  and  bottom  slices. 
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surfaces,  the  principal  stress  difference  is  at  least  greacer  by  a 
factor  of  3.5  than  the  origin  of  line  OA  on  the  top  surface.  In  Fig. 
7->,  the  principal  stress  differences  occurring  along  OB  are  also 
greatest  for  the  bottom.  the  stresses  in  the  middle  were  again  less 
tnan  those  for  top  or  bottom.  in  particular,  the  bottom  stresses  at 
tne  origin  of  OB  were  factors  of  2.P  and  3.5  greater  than  the  top  and 
middle  slices,  respectively.  As  shown  in  Fig.  75,  the  principal 
stress  differences  along  line  OC  are  greatest  for  the  top  slice 
especially  close  to  the  crack  surface.  This  is  understandable  since 
line  OC  is  perpendicular  to  the  intersection  of  the  crack  surface  with 
the  top  surface.  The  stresses  in  the  middle  slice  along  OC  are  still 
less  than  the  top  and  bottom  but  not  as  small  as  in  the  cases  for  lines 
OA  and  OB.  In  fact,  for  the  middle,  stresses  along  OC  are  almost 
identical  to  the  stresses  along  OA  for  the  top  slice.  At  the  trailing 
edge  of  the  crack,  which  corresponds  to  the  origin  of  line  OC,  the 
principal  stress  difference  is  maximum  for  the  top  surface  and  is 
a  factor  of  1.9  greater  than  the  bottom  slice  stresses  and  a  factor 
of  1.7  greater  than  the  middle  slice  stresses. 

A. 5  Mode  I  Opening 

Experimental  data  has  suggested  that  the  crack  front  is 
predominant iv  Mode  I  Opening.  As  shown  in  Fie.  69,  one  of  the 

principal  stress  directions,  is  almost  perpendicular  to  the  fracture 
surface  at  tne  crick  tip.  It  has  been  shown  that  this  as  a  tensile 
stress,  implying  a  Mode  I  loading  is  occurring.  This  would  then 


f 


justify  the  use  of  Mode  I  stress  intensity  factor,  K  in  fracture 
mechanics  calculations  for  the  d.t.  specimen. 
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APPEND  IX  B 

ERROR  ANALYSIS  OF  (V-K)  DATA 

Basically,  two  seeps  were  used  in  calculating  the  error  in  the 
crack  velocity  and  stress  intensity  factor.  The  first  was  the  estimation 
of  the  uncertainties  contained  in  the  equations  for  V  and  K.  These 
uncertainties  were  mostly  due  to  instrumental  errors.  The  second  step 
was  to  use  these  uncertainties  and  employ  the  method  of  propagation 
of  errors  to  calculate  the  final  probable  errors.  It  was  assumed 

in  these  calculations  that  all  the  observed  fluctuations  in  the  data 
were  uncorrelated. 


B . 1  Single  Edge  Notch  Tests 

The  stress  intensity  factor  for  the  s.e.n.  specimen  was  given  by 
Eqn.  (20).  The  uncertainties  for  are  listed  below 

Jp  =  ±  5  N 
Oh/h  =  ±0 . 01 
tw/w  =  ±0.002 
oa  =  ±30 

7f (a/w)/f (a/w)  =  ±0.013 


The  error  in  for  the  s.e.n.  specimen  can  be  written  as 


r;<  -  Ki 

sen 


(?) 2  *  (If) 2  *  (DS) 2  *  (?) 2  if. 


7  ,  ■' 
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For  typical  values  of  a,  and  P,  Eqn.  (54)  yields 


(35) 


The  crack  velocity  for  the  s.e.n.  specimen  was  given  bv  Eqn. 
(19).  The  uncertainties  for  AX  and  At  are  0 ^  =  ±10pm,  =  ±10s. 

The  error  in  V  is  given  by 


°v  * 


(if) 

For  typical  values  in  At,  Eqn.  (56)  yields 


1/2 


(56) 


—  =  *77 

V  ~h° 


(57) 


B . 2 .  Double  Torsion  Tests 

B.2.1  Stress  Intensitv  Factor 


The  stress  intensity  factor  for  the  d.t.  specimen  was  given  ip. 
Eqn.  (7).  The  uncertainties  in  the  parameters  of  the  constant  are 
given  below 


ov/W  =  ±0.0009 
JT/T  =  ±0.025 
cTn/Tn  =  ±0.05 
crWm/Vm  =  ±0.006 
Cv/(l+v)  =  to. 008 
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The  error  in  is  given  by 


+ 


(58) 


/ 


The  uncertainty  in  the  load  measuring  capability  of  the  Instron's  load 
cell  is  approximately  Cp/p  =  ±0.005.  Substituting  these  uncertainties 
inco  Eqn.  (58)  yields 


±4.62 


(59) 


B .  2 . 2  Constant  Deflection  Method  (Decreasing  K^.) 

The  expression  for  the  crack  velocity  during  constant  deflection 
was  given  by  equation  (14).  The  uncertainties  in  the  constant  C9  are 
the  same  as  in  the  constant  except  for  the  uncertainty  in  Young's 
modulus  which  is  given  by  5E/E  =  ±0.05.  The  error  in  is  given  by 


/  CW  \  2 

(OZ)  2 

/  30T  \2  , 

1  2awm  \  “ 

( ov  \2 

1/2 

(t)  + 

It)  + 

(  T  j 

\  wm  / 

(  1+V  ) 

•  (60) 

Substituting  the  appropriate  values  into  Eqn.  (60)  yields 


ar  /C9  *  ±  12.7%  .  (61) 

c2  - 
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The  error  in  after  correcting  for  the  effect  of  plasticity  is 
given  by 


<V  *  %  (fr) 2  +  (f ) 2  +  (r)2  +  25  -  w"6  (6; 

2  9  •? 

For  typical  values,  cry/y  =  0t/t  =  ±  0.002  and  2  £  (P  +  p-)/P“  £  40, 

m  pc 


which  results  in 


oVy/Vy  =  *  13%  .  (63) 

The  major  contribution  to  the  error  in  is  the  uncertainty  in  the 
constant  C2. 

B . 2 . 3  Ccnstant  Load  Method  (Constant  ) 

The  expression  for  the  crack  velocity  during  constant  load  tests 
was  given  by  Eqn.  (18).  The  error  in  after  correcting  for  the 
effect  of  plasticity  is  given  by 

=  [(!r)2  *  iff +  (?)2*  (?) 

The  factor  of  /T”  in  equation  (64)  arises  from  calculating  the  error 
in  the  differences  of  two  quantities  each  with  the  same  error.  The 
ratio  0p/p  for  the  constant  load  method  was  approximately  0.01  and 
using  typical  values  for  the  other  uncertainties  in  Eqn.  (64),  yields 


av  /V 
P  P 


t  18%  . 


(65) 
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VELOCITY  MEASUREMENT  LIMITS 

C .  1  Single  Edge  Notch  Tests 

The  minimum  crack  velocity  that  could  be  measured  using  the 

inscribed  s.e.n.  specimen  and  video  observation  is  limited  onlv  by 

the  patience  of  the  observer.  In  practice,  velocities  as  low  as 

10  ^  m/s  could  be  observed.  The  maximum  crack  velocity  is  limited 

by  the  time  required  for  the  specimen  to  fracture  within  the  duration 

of  one  video  frame.  This  corresponds  to  a  crack  velocity  of  approxi- 

_2 

mately  5  x  10  m/s. 

C . 2  Double  Torsion,  Constant  Deflection  (Decreasing  K.^) 

The  minimum  crack  velocitv  that  could  be  measured  using  the  d.t. 
specimen  depended  on  several  factors.  Among  these  are  the  presence  of 
thermal  fluctuations  in  the  loading  frame,  and  the  stability  and 
precision  of  the  measuring  and  data  acquisition  instruments.  For  the 
constant  deflection  method,  the  limiting  factors  are  the  presence  of 
various  thermal  fluctuations  and  gradients  which  may  affect  the 
measured  load.  The  requirement  to  achieve  testing  temperatures  either 
above  or  below  room  temperature  introduces  additional  thermal 
gradients  in  the  loading  frame  which  may  affect  the  measured  load. 

The  use  of  materials  in  the  loading  frame  with  differing  thermal 
coefficients  of  expansions  wiLl  cause  further  changes  in  the  measured 
load  due  to  ambient  thermal,  fluctuations.  In  this  invest  igation,  toe 
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affects  of  thermal  gradients  were  minimized  by  waiting  a  sufficiently 

long  time  to  achieve  steady  state  conditions  before  testing  a 

specimen.  Thermal  fluctuations  were  minimized  by  air  conditioning 

the  room  and  keeping  the  duration  of  an  average  crack  velocity 

test  from  1/2  to  5  hours.  These  precautions  enabled  crack  velocities 
—  8 

as  low  as  10  m/s  to  measured  with  reasonable  confidence. 

„  .  .  ,  (56) 

Using  precise  temperature  control,  other  investigators  measured 

crack  velocities  as  Low  as  10  ^  m/s  and  Evans  claimed  to  measure 

crack  velocities  down  to  10  ^  m/s  in  glass. 

For  constant  deflection  tests,  upper  crack  velocity  measurements 

were  limited  by  the  response  of  the  load  cell  to  changes  in  load. 

Since  the  load  cell  has  a  high  frequency  response,  crack  velocities  of 

several  hundred  meters  per  second  could  be  measured  reliably.  If, 

however,  a  strip  chart  is  required  to  record  the  load  relaxation,  an 

upper  limit  of  only  several  meters  per  second  could  be  measured. 

C .  3  Double  Torsion,  Constant  Load  (Constant  K j._)_ 

For  the  constant  load  method,  the  limiting  factor  for 
measuring  small  crack  velocities  was  the  stability  of  the.  constant 
load  controller.  Crack  velocities  near  10  m/s  were  measurable  when 
the  load  controller  was  properly  adjusted.  Upper  crack  velocity 
measurement  was  limited  by  the  response  and  speed  of  the  moveable 
crosshead  to  maintain  a  desired  load.  This  placed  a  limit  of 
approximately  5  x  10  m/s  for  the  maximum  velocity  that  could  be 


measured. 


165 


REFERENCES 

1.  H.  K.  Birnbaum,  Environment  -  Sensitive  Fracture  of  Engineering 
Materials ,  Z.  A.  Foroulis,  ed.,  p.  326,  The  Metallurgical 
Society  of  AIME,  Warrendale,  PA,  (1979). 

2.  M.  L.  Grossbeck,  Ph.D.  Thesis,  Univ.  of  Illinois  (1975). 

3.  J.  C.  M.  Li,  R.  A.  Oriani  and  L.  S.  Darken,  Z.  Phys .  Chem.  Neue 
Folge ,  49,  271,  (1966). 

4.  M.  P.  Puls,  Scr.  Me  tall . ,  12:,  595,  (1978). 

5.  M.  P.  Puls,  Atomic  Energy  of  Canada  Limited,  Report  AECL-6302, 
Aug.,  (1978). 

6.  T.  B.  Flanagan,  N.  B.  Mason  and  H.  K.  Birnbaum,  Scr.  Metal  1., 
to  be  published. 

7.  D.  G.  Westlake  and  S.  T.  Ockers,  Metall.  Trans.,  6A,  399,  (1975). 

8.  W.  M.  Albrecht,  M.  W.  Mallet  and  W.  S.  Goode,  J.  Electrochem.  Soc 
105,  219,  (1958). 

9.  H.  K.  Birnbaum,  M.  L.  Grossbeck  and  M.  Amano ,  J.  Less-Common  Met. 
49,  357,  (1976). 

10.  B.  J.  Makenas  and  H.  K.  Birnbaum,  Acta  Metall.,  28^  979,  (1980'. 

11.  M.  L.  Grossbeck  and  H.  K.  Birnbaum,  Acta  Metall.,  2_5  ,  1  33  ,  (  197  7) 

12.  S.  Koike  and  T.  Suzuki,  Acta  Metall.,  to  be  published. 

13.  A.  G.  Evans  and  H.  Johnson,  J.  Mater.  Sci.,  TO,  214,  (1975). 

14.  P.  W.  R.  Beaumont  and  R.  J.  Young,  J.  Mater.  Sci..  1_0,  1334, 
(1975). 

15.  H.  P.  van  Leeuwen,  Corrosion,  31 ,  42,  (1975). 

16.  A.  G.  Evans,  Int.  J.  Fract.,  10,  251,  (1974). 


f 


1  60 


17. 

A. 

G 

.  Evans , 

J.  Mater.  Sci.,  _7>  11  37  ,  (1972). 

18. 

S. 

J 

.  Hudak 

and  R.  P.  Wei,  Metall.  Trans.,  7A,  235,  (1976). 

19. 

R. 

P 

.  Wei,  S 

.  R.  Novak  and  D.  P.  Williams,  Mater.  Res.  Stand 

12 

.  i 

no.  9,  p 

.  25,  (1972). 

20. 

C. 

S 

.  Carter 

,  Corrosion,  2_7>  11,  (1971). 

21. 

H. 

G 

.  Nelson 

,  Metall.  Trans.,  l_k,  621,  (1976). 

22. 

H. 

G 

.  Nelson 

,  D.  P.  Williams  and  J.  E.  Stein,  Metall.  Trans. 

3,  469,  (1972). 


23. 

D. 

P.  Williams  and  A.  G. 

Evans , 

J.  Test.  Eval. 

,  1,  no.  4,  264, 

(1973). 

24. 

D. 

P.  Williams  and  H.  G. 

Nelson, 

Metall.  Trans 

.  ,  1,  63,  (  1970). 

25. 

R. 

A.  Oriani  and  P.  H.  Josephic, 

Acta,  Metall. 

,  25,  979,  (1977) 

26. 

R. 

Dutton,  K.  Nuttall,  M. 

P.  Puls 

and  L.  A.  Simpson,  Metall. 

Trans.,  8A,  1553,  (1977) 

• 

27. 

L. 

A.  Simpson  and  M.  P. 

Puls,  Metall.  Trans., 

10A,  1093,  (1979) 

28. 

L. 

A.  Simpson  and  C.  F. 

Clarke , 

Atomic  Energy 

of  Canada  Limited 

Report  AECL  5815,  Oct., 

(1979). 

29.  J.  0.  Outwater  and  D.  J.  Gerry,  NRL  Interim  Contract  Report, 
Contract  NONR  3219(91) (x),  AD  640848,  Univ.  of  Vermont,  Burlington 
Vt.,  (1966). 

30.  J.  A.  Kies  and  8.  J.  Clark,  Fracture-  1969,  P.  L.  Pratt,  ed., 
Chapman  and  Hall,  Ltd.,  London,  pp.  483-491,  (1969). 

31.  .J.  C.  Pollet  and  S.J.  Burns,  J.  Am.  Ceram.  Soc.,  6_2,  no.  5-6, 
p.  426,  (1979). 

32.  D.  K.  Shettv,  Anil.  V.  Virkar  and  M.  B.  Howard,  J.  .Am.  Ceram.  Soc. 


62,  no.  5-6,  p.  307,  (1979). 


167 


33.  G.  G,  Trantina,  J.  Am.  Ceram.  Soc.,  60,  no.  7-8,  p.  338,  (1977). 

34.  V.  V.  Novozhilov,  Theory  of  Elasticity,  Translation  from  the 
Russian  book  Teoriya  Uprugosi,  Israel  program  for  Scientific 
Translation,  Jerusalem,  p.  232,  (1961). 

35.  0.  W.  Chung  and  N.  S.  Stoloff,  Metall.  Trans.,  9A,  1386,  (1978). 

36.  W.  W.  Gerberich  and  Y.  Katz,  Eng.  Fract.  Mech. ,  1_,  569,  (1969). 

37.  W.  F.  Brown,  Jr.  and  J.  E.  Srawley,  Am.  Soc.  Test.  Mater. 

Spec.  Tech.  Publ.  No.  410,  (1966). 

38.  R.  E.  Peterson,  Stress  Concentration  Factors,  Wiley,  N.  Y.,  (1974). 

39.  S.  J.  Burns,  Private  communications,  (1979). 

40.  E.  R.  Fuller,  Jr. ,  Fracture  Mechanics  Applied  to  Brittle  Materials , 
ASTM  STP  No.  678,  S.  W.  Freiman,  ed. ,  ASTM,  Philadelphia,  PA.  pp . 
3-18,  (1979). 

41.  G.  Matusiewicz,  R.  Booker,  J.  Reiser  and  H.  K.  Birnbaum,  Scr. 
Metall.  ,  8,  1419  ,  (1974) . 

42.  S.  Gahr,  M.  L.  Grossbeck  and  H.  K.  Birnbaum,  Acta  Metall.,  25 . 

125,  (1977). 

43.  H.  Wipt  and  G.  Alefeld,  Phys.  Status  Solidi  (A),  2_3,  1  75  ,  (  197-0. 

44.  G.  Matusiewicz  and  H.  K.  Birnbaum,  J.  Phys.  r.,  7,  no.  11, 

2285,  (1977). 

45.  G.  Schaumann,  J.  Volkl  and  G.  Alefeld,  Phvs.  Status  Solidi  (A), 

42,  401,  (1970). 


46. 

D.  I. 

Bolef , 

J .  Appl .  Phys  .  ,  32 , 

100,  (1961). 

•4  /  . 

M.  S. 

Rashid 

and  T.  E.  Scott,  J. 

Less-Common  Met 

168 


48.  C.  C.  Chen  and  R.  J.  Arsenault,  Hydrogen  in  Metals,  I.  M. 
Bernstein  and  A.  W.  Thompson,  ed.,  p.  393,  Metals  Park,  Ohio, 
American  Society  for  Metals  (1974). 

49.  I.  A.  Pereturina  and  V.  A.  Pavlov,  Phys.  Met.  Metallogr.,  41 , 

No.  4,  p.  146,  (1976). 

50.  S.  Gahr,  B.  J.  Makenas  and  H.  K..  Birnbaum,  Acta  Metall.  ,  28, 
1207,  (1980). 

51.  A.  W.  Hendry,  Elements  of  Experimental  Stress  Analysis,  S.  I. 
Edition,  Pergamon  Press,  Oxford,  p.  97,  (1977). 

52.  A.  J.  Durelli  and  W.  F.  Riley,  Introduction  to  Photomechanics, 
Prentice-Hall,  Inc.,  Englewood  Cliffs,  N.J.  ,  p.  64,  (1965). 

53.  S.  P.  Timoshenko  and  J.  N.  Goodier,  Theory  of  Elasticity,  third 
Edition,  McGraw  Hill,  N.Y.,  p.  167,  (1970). 

54.  H.  T.  Jessop  and  F.  C.  Harris,  Photoelasticity,  Principles  and 
Methods,  Dover,  N.  Y.  ,  p.  73,  (1960). 

55.  P.  R.  Bevington,  Data  Reduction  and  Error  Analysis  for  the 
Physical  Sciences,  McGraw  Hill,  N.Y.  ,  p.  56,  (1969). 

56.  B.  J.  Pletka,  E.  R.  Fuller,  Jr.  .  j.  G.  Koepke,  Fracture 

Mechanics  Applied  to  Brittle  ,  ASTM  STP.  No.  6  78,  S.  W. 


Freiman,  ed.,  ASTM,  Philadelphia,  PA,  pp.  19-37  (1979). 


169 


VITA 

Barry  Scott  Hindin  was  born  on  January  19,  1953  in  Chicago, 
Illinois  and  resided  there  until  graduating  from  Stephen  Mather  High 
School  in  1970.  As  a  junior  high  school  student  in  1969,  he  entered 
the  Chicago  Public  School  Science  Fair  contest  and  his  paper  entitled 
"The  Effects  of  Electromagnetic  Fields  on  the  Acceleration  of  Steel 
Ball  Bearings"  won  first  place  in  the  physics  division.  He  enrolled 
at  the  University  of  Wisconsin  at  Madison  in  1970  and  as  a  physics 
major  established  and  became  president  of  the  undergraduate  Physics 
Club  in  1973  and  graduated  with  a  Bachelor  of  Science  Degree  in  1974. 
Since  1974  he  has  been  a  graduate  research  assistant  in  the  Department 
of  Metallurgy  and  Mining  Engineering  at  the  University  of  Illinois  at 
Urbana-Champaign.  In  August  1979,  he  married  Emily  Louise  Sudman. 


{ 


DOCUMENT  CONTROL  DATA  -  R  &  D  j 

'Security  c  fas  at  fieation  of  title,  body  of  sbstrect  end  indewing  annotation  muii  6a  entered  when  the  overall  report  Is  classified) 

\  ORi  G .  n  a  ting  *c  TIVI  Tv  (Corporal#  aufftorj 

University  of  Illinois  at  Urbana  Champaign 

2e.  REPORT  security  classification 

Unclassified 

16.  GROUP 

1  Tl  TLfe 

Fracture  Kinetics  of  Hydrogen  Embrittled  Niobium 

•«  Descriptive  NOT€J  (7\pe  of  report  end  inclusive  detea) 

Technical  Report 

S  a Vj  ThORUl  (First  name,  middle  initial,  laet  name) 

Barry  s  Hindin  and  H  K  Birnbaum 

1  6  Rt  PORT  OATE 

|  March  1981 

7*.  TOTAL.  NO  O  F  PAGES  7b.  NO  OF  REFS 

168  Sfi 

*«<  contract  or  grant  no 

USN  00014-75-01012 

o.  *» roj e:t  no 

c. 

d. 

9a.  ORIGINATOR'S  REPORT  N  L>  M  6  £  R(S) 

96.  OTHER  REPORT  NOI3I  or.Ntr  numfavrs  that  may  b*  assigned 

this  report) 

IC  DISTRIBUTION  JTATEMENT 

This  document  is  unclassified.  Reproduction  and  distribution  for  any  purpose 
of  the  US  Government  is  permitted. 

M  JjPOUMENTARV  NOTES 

»2.  sponsoring  military  activity 

Office  of  Naval  Research 

3  ABSTRACT 


See  attached 


Unclassified 

Security  Cla^si  ficjnon 


DD  ?::.,1473 


